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ABSTRACT 

We present the chemistry, temperature, and dynamical state of a sample of 193 dense cores or core 
candidates in the Perseus Molecular cloud and compare the properties of cores associated with young 
stars and clusters with those which are not. The combination of our NH3 and CCS observations 
with previous millimeter, sub-millimeter, and Spitzer data available for this cloud enable us both 
to determine core properties precisely and to accurately classify cores as starless or protostellar. 
The properties of cores in different cluster environments and before-and-after star formation provide 
important constraints on simulations of star-formation, particularly under the paradigm that the 
essence of star formation is set by the turbulent formation of prestellar cores. We separate the influence 
of stellar content from that of cluster environment and find that cores within clusters have (1) higher 
kinetic temperatures (12.9 K vs. 10.8 K) and (2) lower fractional abundances of CCS (0.6xl0~^ 
vs 2.0x10"^) and NH3 (1.2xl0~® vs 2.9x10^^). Cores associated with protostars have (1) shghtly 
higher kinetic temperatures (11.9 K vs. 10.6 K) (2) higher NH3 excitation temperatures (7.4 K vs. 
6.1 K), (3) are at higher column density (1.2x10^^ cm"^ vs. 0.6x10^^ cm~^), have (4) shghtly more 
non-thermal/turbulent NH3 linewidths (0.14 km/s vs. 0.11 km/s FWHM), have (5) higher masses 
(1.5 M0 vs. 1.0 M0) and have (6) lower fractional abundance of CCS (1.4x10"^ vs. 2.4x10"^). 
All values are medians. We find that neither cluster environment nor protostellar content makes a 
significant difference to the dynamical state of cores as estimated by the virial parameter ~ most cores 
in each category are gravitationally bound. Only the high precision of our measurements and the size 
of our sample makes such distinctions possible. Overall, cluster environment and protostellar content 
have a smaller influence on the properties of the cores than is typically assumed, and the variation 
within categories is larger than the differences between categories. 
Subject headings: ISM:clouds — ISM: molecules — radio lines:ISM 



L INTRODUCTION 

Stars form in dense cores, providing ample motiva- 
tion for understanding these enigmatic objects. Several 
authors have claimed that the stellar initial mass func- 
tion (IMF) comes directly from the core-mass function 
(CMF) with a Saltpeter power-law slope at high masses 



population are required to constrain theoretical models. 

The deflnition of a dense core differs depending on the 
method employed, sensitivity achieved, and physical res- 
olution reached. The platonic ideal is a small (0.1 pc), 
cold (~ lOK), dense (> 10^ cm~^), quiescent (smal l non- 



thermal linewidth) object ([Myers fc Benson 1983 [Ben 



(dN(m)/d log{m) ~ m-i'^S) and a turnover at some |son fc Myers| [l989| [Goodman et al.| |1998p witFiltrong 



characteristic mass, although the valu e of this mass is 
in dispute (compare Motte et al. ( 1998 1 with Alves et al. 
([2007 )). This correspondence between the mass tunc- 
tions suggests two important corollaries: that star forma- 
tion is a simple process which transforms core mass into 
stellar material in a robust and largely invariant fash- 
ion; and that a theory of molecular cloud formation and 
evolution which provides the correct core-ma ss function 
(e.g. the turbulent fragmentation model of Padoan & 
Nordlund (2002)) has solved a critical part of the star 
formation problem. However, Swift & Williams (20081 
model several diverse evolutionary schemes which map 
a CMF onto an IMF and show that the present data on 
mass functions is insufficient to distinguish between their 
different schemes. Other properties of the dense core 
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density contrast from its surroundings and no significant 
substructure, allowing a si mple radial power-law descrip - 
tion of their structure (e.g. Ward-Thompson et al.|1994 l. 
In addition, if we are searching tor the precursor to stars 
we want to exclude transient objects; a dense core should 
be close to gravitationally bound. However, gravitation- 
ally unbound cores m ay be compressed by turbulence 
and thus still collapse ( Gomez et al.|2007 1 . It is also not 
clear whether we should exclude a dense core which con- 
tains a young star. We know that such an object is the 
formation site of a star, but the protostar may have had 
significant impact upon the core, rendering it useless for 
studying the initial conditions of star formation. 

Many detailed studies have been carried out on indi- 
vidual dense cores (see Bergin & Tafalla (20071 for a 
review) . Much can also be gained from cloud-wide stud- 
ies of cores, for though they must sacrifice some level of 
detail, such studies are able to average over certain in- 
evitable biases (such as lacking velocity information on 
the plane of the sky and spatial information along the line 
of sight), provide good statistics of the population as a 
whole, and allow us to assess the infiuence of environmen- 
tal variations. One particularly powerful combination is 
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that of NH3 observations with either continuum emis- 
sion (typically millimeter or sub-millimeter) or extinc- 
tion maps which trace the dust within cores. Sufficiently 
high signal-to-noise ammonia spectra provide a measure- 
ment of the temperature and the total linewidth of an ob- 
ject, so we are able to say whether an object is cold and 
quiescent. The critical density of ammonia is roughly 



et al. 20051 surveys of Perseus. Follow-up studies have 
examined the distinction betwe en mass distributions of 
protostellar and starless cores dHatchell fc Fuller||2008 



Enoch et al.|2008| |Kirk et al.|20b6p , so we will touch but 
lightly on this question, 'i'he stability or dynamical state 



of these cores have also been previously studied ( Kirk 



et al. 2007 



10 cm ^ (Ho & Townes 19831, so objects detected in for a more detailed analysis. 



Enoch et al. 2008 1 , but our data sets allow 



NH3 are typically at or above this density and thus more 
dense than the bulk of the molecular cloud. Ammonia 
also provides one piece of the stability puzzle, namely 
how much internal non-thermal velocity dispersion other 
forces must overcome in order to collapse the core. Mil- 
limeter maps can provide spatial extent, disambiguation 
from the background, and masses (with the aid of mea- 
sured or assumed temperatures). This provides a first 
order estimate of whether a core is bound, though mag- 
netic support and external pressure still elude us. Stellar 
content can be assessed through direct infrared observa- 
tions of point sources with an instrument like Spitzer. 

We present such a synthesis of observations for the 
molecular cloud in Perseus by combining the results of 
an ammonia sur vey of core-candidates w ith the Green 



The distinction between clustered and isolated cores 



Bank Telescope (Rosolowsky et al. 2008|) with a Bolo 



cam survey of the continuum at 1 mm ('Enoch et al. 2006 1 



and a catalog of young stellar objects from the C ores-to 
Disks Spitzer Legacy project ( Evans et al.||2003" l. These 
three data sets are all cloud- wide and spatially complete. 
The NH3 survey covers every dense core and dense core 
candidate within the cloud, drawn principally from the 
Bolocam survey with addi tional candidates drawn from 
additional COMPLETE^ (|Ridge et alj[2006| data as de- 
scribed in [Rosolowsky et 1 2008^ 

Our analysis includes fewer cores (of the 193 poi ntings, 
122 are Bo locam cores) than the seminal work of Jijina 
et al. ( 1999) (264 cores), but is complete for Perseus down 
to the mass limit of the Bolocam survey, and is more 
uniform. This uniformity comes both from using a con- 
sistent observational setup and from the fact that we are 
studying one cloud at a fixed distance 



The study of the Pipe Nebula in Lada et al. ( 2008 1 is 



a comparable cloud-wide study based on extinction map 
made from 2MASS, C^^ O observations from t he Arizona 
Radio Observatory 12m dMuench et"Sl|2007|, and NH3 
observations from the GBT ( [Rathborne et al.|2008| . Be- 
cause Perseus has more massive cores and is at a better 
declination for observing from the GBT, we are able to 
detect NH3 from every Bolocam core. This allows us to 
measure the temperature of each core. The Pipe also has 
extremely limited star formation, and may represent an 
early stage of star formation, while Perseus provides ex- 
amples of both isolated and clustered star formation as 
well as relatively pristine starless objects. 

This paper uses a statistical approach to separate and 
measure the influence of both stellar content and clus- 
ter environment on dense cores. Our main focus will 
be on the chemistry, temperature, and dynamical states 
of the cores. The mass distributions of cores through- 
out Perseus have already been explored in detail by a 
number of pa pers using two the rmal emission surveys: 
the Bolocam ( [Enoch et al]|2006| ) and SCUBA ( [Hatctiell 



^ coordinated Molecular Probe Line Extinction Thermal Emis- 
sion Survey of Star Forming Regions 



was made in [Jijina et al. (19991, but has not since re- 
ceived sufficient attention; our cloud-wide survey is the 
perfect data-set for such a study, as Perseus features two 
main clusters as well as star formation i n isolated regions- 



Recen t models of star formation (e.g. Nakamura & Li 



2007) which include the role of feedback clearly estab 
lish the need to assess the influence of environment upon 
the star formation process, an influence which may man- 
ifest itself in the core properties. If cores within clusters 
are substantially different from isolated cores, this may 
explain why the CMF appears to be different in differ- 
ent regions. It is also possible that the role of feedback 
means that star formation within clusters is a different 
process than star formation in isolated regions and that 
the notion of a simple function which maps the CMF to 
the IMF must at least be modified to allow for different 
environments. 

We describe our datasets in Sj2] and describe our sta- 
tistical tests in ^ In Q we present in full detail the 
distributions of the following variables: kinetic tempera- 
ture, excitation temperature, column density, fractional 
abundance of NH3 and CCS, linewidths (total and non- 
thermal), mass, and stability (the virial parameter, a). 
Along the way, we compare with results obtained for 
Perseus in previous studies, and present a comparison 
with studies in other regions in ^ Distribut ions of all 
variables by category are provided in Table |2] while a 
summary of which variables are significantly influenced 
by cluster environment or protostellar content is provided 
in Table [S] 

2. DESCRIPTION OF DATASETS 
2.1. Ammonia Survey 

The ammonia data in this paper are taken fro m the 
spectral atlas described in Rosolowsky et al. (2008 1. The 
survey consists of 193 spectra taken with the Green Bank 
Telescope, pointed towards a selection of cores and core 
candidates in the Perseus molecular cloud. The targets 
were drawn from a merged list designed for maximum 
completeness based on several surveys. In addition to 
cataloged sources of millimeter continuum emission, the 
survey targeted marginal millimeter sources and compact 
regions of extinction. The GBT spectrometer was config- 
ured to produce high resolution (0.025 km s""'^) spectra of 
NH3 (1,1) and (2,2) emission as well as CCS (2i Iq) 
and CC^^S (2i -> Iq). The beam size of the GBT is 
31"at these frequencies which projects to a size of 0.04 
pc at the distance of Perseus (250 pc, see Enoch et al. 
2006 for a full discussion of the distance). 

Fmlowing reduction, a uniform slab radiative transfer 
solution was fit to the full ammonia complex to determine 
the values of the kinetic temperature (Ty^in), the line-of- 
sight velocity (wisr), the velocity dispersion of the tracer 
((Tnhs), the radiative excitation temperature (Tex) and 
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the full line optical depth (r). Using a combination of the 
derived parameters, the column densities of both NH3 
and CCS were also derived. Since the model is necessar- 
ily simple, complicated sources will be represented by av- 
erage properties. However, the model produces good fits 
for the wide variety of sources found in the survey data. 
In cases where the spectra were clearly separable into 
two distinct components such fitting was performed, but 
these results were fiagged and excluded from our anal- 
ysis of variables which rely upon an integrated column 
measurement (fractional abundances, dust masses, and 
the virial parameter). 

2.2. Bolocam Survey 

For a fiducial sample of dense cores, we used the mil- 
limeter c ontinuum observat ions of dust emission pre- 
sented by Enoch et al. (2006 1. That work used Bolocam, 
a bolometer array on the Caltech Submillimeter Observa- 
tory, to map the high extinction portion {Ay > 2) of the 
Perseus molecular cloud at A « 1.1 mm. Using typical 
assumptions about dust emission, the Bolocam survey is 
complete down to ~ 0.2 Mq at the assumed distance of 
Perseus. However, owing to variations in the degree to 
which sources are recovered under the Bolocam mapping 
strategy, the mass sensitivity is governed by the source 
size. In addition, large scale structure is filtered out by 
the adopted observation strategy, though this filtering 
tends to highlight the dense cores desired for study. 

Dense cores are identified by convolving the final map 
by an optimal filter that selects point sources. From a 
list of significant peaks in the filtered map, only sources 
with a well defined centroid are considered dense cores. 
Fluxes are measured in apertures around each object and 
major/minor axis sizes are from two-dimensional Gaus- 
sian fits to the emission profile. A total of 122 dense 
cores are so identified and cataloged. These cores were 
t he p rimary catalog adopted for the spectral line survey 
(j |2.1[ ) owing to the complete spatial coverage provided 
by the Bolocam survey and the similarity in resolution 
to the GBT (also 31"). The Bolocam survey likely iden- 
tifies every region harbo ring dense cores, though higher 
resolution observations ( [Kirk et al. 2006 1 suggest that 
some of the cores may have significant substructure. 

Although there is evidence that in certain high-density 
and high-radiation environments NH3 does a poor job 
of tracing some dust cores (see the study of Friesen et 
|al I (|2009| in Ophiuchus), maps of regions ot Perseus 
made with GBT NH3 observations similar to those in 
this pointed survey show an excellent correspondence 
with dust ma ps such as the BOLOCAM map used in 
this study (Pineda et al.||2009[ ). 



2.3. YSO and Cluster identification 

Classifying cores as starless or protostellar is a difficult 
problem. In the case of Perseus, several previous stud- 
ies have examined this issue from diffe rent perspectives 
( Hatchell et al.|2007 Enoch et al.|2008 ) , so our approach 
is to synthesize these efforts with some of our own meth- 
ods. The areas of disagreement between these methods 
hint at the complexity of this problem which is exac- 
erbated by differing definitions as well as observational 
uncertainty. Full details of this synthesis are presented 
in Appendix A, but we take the abundantly cautious ap- 
proach of tagging each core as either (1) Starless, (2) 



TABLE 1 

Number of cores in each category 

Starless Protostellar 
certain(uncortain) certain(uncortain) 



Clustered 
Isolated 



23 (7) 
102 (4) 



16 (3) 
38 (6) 



Subset of cores with NH3 detections 



Clustered 
Isolated 



19 (7) 
79 (4) 



15 (3) 
35 (6) 



Probably Starless, (3) Probably Protostellar, or (4) Pro- 
tostellar, with the uncertain classifications resulting from 
disagreements between the various methods used. This 
classification allows us to verify that cores with uncertain 
classification do not influence our results. Our final list of 
starless/protostellar for these NH3 pointings is presented 
in Table [5] in Appendix A, with comments. 

We classified each core as isolated or clustered based 
on the definitions of the cluster boundaries of IC 348 
and NGC 1333 in lj0rgensen et al.J ( |2006| , identified 
as high concentrations of YSO candidates from Spitzer 
data. For NGC1333 this means from 03:28:00 to 03:30:00 
in RA and -^31:06:00 to -1-31:30:00 in Dec, while IC 348 
goes from 03:43:12 to 03:46:00 in RA and -1-31:48:00 to 
+32:24:00 in Dec. For NGC 1333 there is a clear drop 
of in the density of objects beyond this boundary, while 
in IC 348 the cut-off at the western edge is slightly more 
ambiguous, but different definitions of the clusters would 
change only a small number of objects. 

The number of cores within each category is presented 
in Table [T] These numbers represent the upper limits on 
how many objects may be considered in our analyses. For 
each variable, we use only objects for which that variable 
is well measured. In particular, some objects lack either 
significant NH3 or millimeter emission and are unlikely 
to be genuine cores. 

If we restrict ourselves to objects with a NH3 detec- 
tion, the ratio of starless to protostellar cores is higher in 
isolated regions (Nstariess/Nprotostciiar = 2.0) than in clus- 
tered regions (Nstaricss/Nproto steiiar = 1.4). Both t hesc ra- 
tios are higher than found in |Enoch et al.| (|2008| , where 
(Nstarioss/Nombodded = 1-0) for Pcrscus. This is partly due 
to our different classification of which cores are starless as 
explained in Appendix A, but is mostly due to the extra 
71 positions we target which are not in the Bolocam sam- 
ple and are largely starless. The diverse datasets from 
which these additional positions were drawn prohibit us 
from making a strong inference from these ratios. 



3. STATISTICAL METHODS 

3.1. Distributions of Variables 

The following variables are derived as part of the fit o f 
the NH3 and CCS spectra in Rosol owsky et al.| ( |2008[ ): 
Tkin, Tex, NH3 and CCS column density, NH3 andCCS 
total linewidths. In this paper, we use Tkin from the 
NH3 line and Bolocam fluxes derived in lEnoch et al.| 



( 2006 1 in order to derive total column densi ties (j |4.3[ ) arid 
fractional abundances of NH3 and CCS (! |4.4[ ). We use 
Tkin and the total linewidths to calculate non-thermal 
hnewidth of both NH3 and CCS (p!5|. Finally, we use 
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Tkin to improv e the estimated ma sses (j |4.6[ ) of Bolocam 
cores derived in Enoch et al. (2006 1 and from this and the 
n on-t hermal linewidths we estimate the virial parameter 
(j |4.6[ ) for each core. The derivation of new quantities is 
presented in the sections referenced above. 

Table [2] presents a summary of the distributions of all 
our variables, broken down by classification. Values are 
presented for the mean, standard deviation, median, and 
25% and 75% quartiles, as many variables are not well 
described as a normal distribution. As we analyze only 
cores where a given variable can be well determined, we 
suffer some bias due to non-detections of small values. 
This is possibly a concern in Ty^i^ and is definitely a con- 
cern for the fractional abundance of CCS (X-CCS), mass, 
and the virial parameter, a. We discuss these biases 
within the text for each variable. 

Overall, the variation within categories is generally 
larger than the differences between categories and only 
the high precision of our measurements and the size of 
our sample makes it possible to draw distinctions be- 
tween categories. Table [3] highlights the distributions 
where statistically significant differences exist, as deter- 
mined by an Analysis of Variance test, described below. 

3.2. Analysis of Variance Statistical Test 

In order to quantify the respective influence of stellar 
content and cluster environment, we performed Analysis 
of Variance tests (ANOVA ) on a numbe r of key param- 
eters. ANOVA tests (see Scheffe 1959 for a complete 
reference) are useful for assessmg the influence of cate- 
gorical or qualitative variables on other variables. In this 
case, we categorize each core as starless/protostellar and 
clustered/isolated and test the hypothesis that these cat- 
egories have a significant influence on the measured quan- 
titative variables. In particular, the Two- Way ANOVA 
framework allows us to construct a statistical model to 
test for the "interaction effect" by comparing a model 
containing an interaction term: 



against a model without this term. 



y^jk ^ fi + ai + (3j + 



Hjk, 



(1) 



(2) 



via the F-test. In these models we have two variables, 
a at i levels and /3 at j levels. We test the hypothe- 
sis that each individual core's property (Y^) is due to 
some mean value (/x), some contribution from protostel- 
lar content (a), some from environment (/3), some from 
the interaction effect (a/3), and some from random error 

A simple hypothetical example of the interaction ef- 
fect would be the following: starless cores are 2 K colder 
than protostellar cores, and isolated cores are 1 K colder 
than clustered cores, but starless AND isolated cores 
are 5 K colder than protostellar cores in clusters - i.e. 
the influence of the two variables is not simply addi- 
tive/linear. In our study it turned out that in all cases 
the interaction effect was statistically insignificant, which 
allows us to make direct statements about the influence of 
the individual categories (starless/protostellar and clus- 
tered/isolated). For the starless/protostellar distinction, 
we treated uncertain classifications as distinct categori- 
cal variables. The separate classification of uncertain ob- 
jects essentially eliminated these objects from our anal- 



ysis, but allowed us to check for any unusual properties 
these uncertain objects might have had. 

As with many traditional statistical tests, ANOVA as- 
sumes that the underlying variables are normally dis- 
tributed, but this is generally not the case with our mea- 
sured variables. However, in many cases the variables 
can be transformed with a logarithmic transformation to 
be approximately normal, and this allows us to inter- 
pret the statistical significance of our findings within the 
ANOVA framework. Differences between means of loga- 
rithmically transformed variables are expressed as ratios 
of the median values because the transformation is mono- 
tonic (preserves the median) and the transformed distri- 
bution is roughly normal so the median approximates the 
mean. The differences between the category variables are 
generally small compared to the spread within individual 
populations. 

To derive realistic confidence interva ls we used Tuk ey's 
Honest Significant Differences (HSD) ( |Yandell||1997 pp. 
97-98), which account for the fact that we were approach- 
ing the data set without specific comparisons to test, but 
were rather testing the significance of all possible pair- 
ings. This is the most conservative approach to setting 
confidence intervals. 

Our ANOVA procedure makes the assumption that the 
variance within groups are the same. We tes ted this as- 
sumption using Levene's test (Levene 19601, and cases 
where significant non-uniformity where encountered are 
mentioned in the text and summary table. These signif- 
icance levels and confidence intervals should be treated 
with greater suspicion, but the ANOVA tests are rel- 



atively insensitive to non-uniform variance (see [Levene 
19601, so the general sense of the result is probably un- 
attected. 

We performed these tasks in the R statistical package. 
The results are summarized in Table [3] and will be dis- 
cussed in more detail for each variable in the following 
sections. 

3.3. Sensitivity to Random and Systematic 
Temperature Errors 

Most of the derived variables we study herein rely in 
some way on using a temperature (Tkin) derived from 
modeling the NH3 spectra to convert observed quantities 
into physical units. This dependence on temperature can 
be divided into two cases: 

1) The distributions which use only the NH3 and CCS 
data (Tox, Tkin, and both non-thermal linewidths) re- 
quire that our determination of temperature accurately 
measures the temperature of the bulk (H2) gas. The 
most serious error possible is a systematic one. 

2) The distributions which also rely on the Bolocam 
maps (N(H2), Mdust, a, X(CCS) and X(NH3)) require 
that Tkin from NH3 adequately approximates Tjust- This 
is an estimate and may be biased in some way. We show 
in Appendix B that the best estimate for this discrepancy 
has a small influence on our results. Nonetheless, for 
this half of the analysis the probabilities quoted are valid 
only under the assumption that Tkin for NH3 adequately 
approximates Tdust- 

As the two simplest quantities to assess, we perform 
experiments on the non-thermal linewidth of NH3 (case 
1) and the total column density (case 2) in Appendix 
B. These tests show that our results are very robust 
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TABLE 2 

Summary of Variable Distributions by Classification 



Variable 


Classification'' 


Mean 


a 


Median 


25% quartile 


75% quartile 


Tkin (K) 

Kinetic temperature 


Protostellar 
Starless 
Cluster 
Isolated 


12.67(12.51) 
11.43(11.75) 

13.35 

11.51 


2.4(2.3) 
2.7(3.1) 

3.2 

2.5 


11.85(11.83) 
10.55(10.67) 

12.92 

10.79 


11.32(11.01) 
9.99(10.02) 

11.15 

10.08 


14.05(13.49) 
11.66(12.41) 

14.40 

11.96 


Tex (K) 

Excitation temperature 


Protostellar 
Starless 
Cluster 
Isolated 


7.25(7.12) 
5.99(6.13) 

6.82 

6.48 


1.2(1.3) 
0.8(0.9) 

1.3 

1.2 


7.35(7.26) 
6.09(6.13) 

6.57 

6.41 


6.46(6.19) 
5.58(5.74) 

6.03 

5.83 


8.32(8.08) 
6.57(6.60) 

7.82 

7.35 


N(H2) (lxl022 cm-2) 
Total column density 


Protostellar 
Starless 
Cluster 
Isolated 


2.22(2.08) 
0.66(0.73) 

1.13 

1.23 


2.7(2.5) 
0.5(0.6) 

1.0 

1.9 


1.18(1.18) 
0.59(0.63) 

0.93 

0.70 


0.70(0.71) 
0.28(0.28) 

0.52 

0.44 


2.94(2.59) 
0.85(0.93) 

1.42 

1.24 


X-NH3 (1x10-8) 
Fractional abundance of NH3 


Protostellar 
Starless 
Cluster 
Isolated 


5.39(5.15) 
8.17(7.61) 

1.62 

8.57 


13.5(12.5) 
32.3(31.0) 

1.4 

30.8 


3.39(3.39) 
2.31(1.87) 

1.18 

2.89 


1.06(1.09) 
0.63(0.52) 

0.52 

0.86 


5.00(5.04) 
4.40(4.16) 

2.59 

5.19 


X-CCS (1x10-8) 
Fractional abundance of CCS 


Protostellar 
Starless 
Cluster 
Isolated 


1.48(1.83) 
3.24(3.11) 

0.86 

2.87 


0.9(2.1) 
2.5(2.5) 

0.4 

2.5 


1.43(1.43) 
2.37(2.32) 

0.60 

1.96 


0.79(0.79) 
1.70(1.58) 

0.55 

1.43 


2.06(2.06) 
4.65(4.32) 

1.50 

3.24 


f^non-thcrmal (NH3) (km/s) 

Non-thermal linewidth 


Protostellar 

Starless 
Cluster 
Isolated 


0.142(0.138) 
0.117(0.122) 

0.138 

0.124 


0.04(0.04) 
0.04(0.04) 

0.04 

0.04 


0.135(0.134) 
0.114(0.117) 

0.135 

0.122 


0.113(0.113) 
0.097(0.099) 

0.114 

0.103 


0.176(0.174) 
0.137(0.148) 

0.169 

0.145 


CTnon-thermal (CCS) (km/s) 

Non-thermal linewidth 


Protostellar 
Starless 

V_^1L1.0 LCI 

Isolated 


0.196(0.195) 
0.165(0.172) 

0.135 

0.186 


0.09(0.08) 
0.08(0.08) 

0.06 

0.08 


0.181(0.181) 
0.151(0.155) 

0.129 

0.165 


0.153(0.153) 
0.120(0.120) 

0.129 

0.144 


0.226(0.226) 
0.180(0.210) 

0.155 

0.217 


Mass (M0) 

Mdust 


Protostellar 
Starless 
Cluster 
Isolated 


2.37(2.44) 
1.41(1.35) 

1.57 

1.91 


2.8(3.0) 
0.9(0.9) 

1.3 

2.3 


1.50(1.46) 
1.02(1.02) 

1.14 

1.20 


0.76(0.76) 
0.82(0.82) 

0.76 

0.81 


2.80(3.45) 
1.90(1.81) 

1.92 

2.30 


Virial parameter: a 
^dynamical /Mdust 


Protostellar 

Starless 
Cluster 
Isolated 


1.36(1.31) 
1.37(1.35) 

1.47 

1.28 


0.8(0.8) 
0.5(0.5) 

0.8 

0.6 


0.99(0.99) 
1.39(1.39) 

1.37 

1.24 


0.80(0.75) 
1.12(1.04) 

0.82 

0.84 


1.97(1.97) 
1.64(1.58) 

2.06 

1.57 



^ For the protostellar/starless classifications, the results are displayed as: certain(including uncertain). 



against the influence of random errors on the tempera- 
ture. This is due to the high-quality of our spectra, which 
allow a very precise estimate of the temperature. The 
non-thermal linewidth result is robust against system- 
atic errors in the temperature determination. System- 
atic errors in the temperature determination can be con- 
structed which eliminate the distinctions we see in quan- 
tities which rely on the dust emission map, particularly 
if these systematic differences work in different directions 
for the two populations (i.e. protostellar/starless). How- 
ever, under our best estimate of the possible systematic 
error our results are aifecited by only a small amount. 

We use these experiments to infer the influence of ran- 
dom and systematic temperature errors on our other vari- 
ables. 

The non-thermal linewidths of CCS will behave much 
like the non-thermal linewidths of NH3; they will be in- 
fluenced only slightly by either sort of temperature error. 

Errors in total column density will propagate into the 
fractional abundances of NH3 and CCS mostly in the 



obvious way (increasing the total column density of ma- 
terial decreases the fractional abundance). A bad flt to 
the NH3 could trade off some NH3 column for some Ttin, 
but as our spectra are generally of high quality and our 
fits are generally quite good, we do not think this is likely 
to be a problem. Masses (which come directly from the 
dust map) will obviously be affected in the same was a 
column density. The virial parameter (a = Afdyn/Mdust) 
has a double dependence on Tkin- However, since our ex- 
periments in Appendix B show non-thermal linewidths to 
be influenced by only a small amount by reasonable er- 
rors on Tkin; the doiniiiaiit effec;t will be in Mjust- Thus 
if systematic errors are present which cause N(H2) to be 
overestimated for a given object, a will underestimated, 
and vice-versa. 

4. RESULTS 

4. 1 . Temperature 
4.1.1. Kinetic Temperatures: Tkin 
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TABLE 3 
Summary of ANOVA Results 



Variable 




Cluster Signifieant? 


Protostar Significant? 


Statement 


Oomments 


# of Cores 


Tkin 


log 


Yes, P = 6x10-5 


0.01 


Med(TA') Starred:Starless = 1.11 
95% 01(1.01,1.23) 

Mecl(A.-UCo) Uiusteriisolated = i.lD 
95% 01(1.08,1.24) 


Long-tailed after 
log transform 


128 


Tex 




No 


Yes, P = IxlO"'^ 


Tex Starred - Starless= 1.26 K 
95% 01(0.73,1.79) 


Normal, unequal 
variance'' 


118 


N(H2) 


log 


No 


Yes, P = 3x10^10 


Med(N(H2)) Starrcd:Starless = 2.82 
95% 01(1.92,4.16) 


Normal after 
log transform 


181 


X-NH3 


log 


Yes, P = 0.02 


No 


Med(X-NH3) Cluster;Isolated = 0.56 
95% 01(0.35,0.91) 


Skewed after 
log transform 


178 


X-CCS 


log 


Yes, P = 8e-8 


Yes, P = 2x10-5 


Med(X-OCS) Starred:Starless = 0.37 
95% 01(0.22,0.62) 

Med(X-OCS) Cluster:Isolated = 0.32 
95% 01(0.21,0.48) 


Uses Icr limits 
as data points. 
Similar results 
without limits 


181 


Vnt (NH3) 




No 


Yes, P = 4x10-3 


Vnon- thermal Starred - Starless = 

0.025 km/s 

95% 01(0.003,0.048) 


Normal 


83 


Vnx (CCS) 




No 


No 


No significa-nt influence 


Still long-tailed, 
no COS in proto- 
stellar -|- cluster 


59 


Mass 


log 


No 


Yes, P = 0.05 


Med(Mass) Starred;Starless = 1.32 
95% 01(1.00,1.76) 


Normal after 
log, unequal 
variance'^ 


123 


a 


log 


No 


No 


No significant influence 


Normal after log 
transform 


69 



^ The transformation applied to make the distribution more normal. All logs are log base 10. Using Lcvene's test we find the cxeitation temperatures 
of protostellar and starless eores to have unequal variance at the P — 0.02 level. Using Levenc's test wc find the log(Masses) of protostellar and starless 
cores to have unequal variance at the P — 0.01 level. 



The kinetic temperature of cores were fit from the (1,1) 
and (2,2) transitions of NH3 and depend on the relative 
strengths of these two transitions. They are thus deter- 
mined only for objects with a detected (2,2) line. Our 
ability to detect this line depends on the NH3 column 
and temperature of the core. Thus, we are generally in- 
sensitive to very low kinetic temperatures, and at low 
temperatures we are biased against being able to solve 
for the temperature of diffuse or low-density material. 
This may not be a significant problem for our survey, as 
we illustrate in Figure [l] which shows T^in versus the 
integrated intensity of the (1,1) line. Good solutions for 
Tkin are shown as plus symbols and 3 a upper limits on 



Tkin (derived in Rosolowsky et al. |2008[ using the noise 
in each individual spectrum) are shown for the rest of the 
objects. These upper limits are consistent with our sur- 
vey not being significantly biased against detecting cores 
with Tkin below 9 K. That is, the (1,1) is always quite 
weak for objects without a (2,2) detection. We cannot 
definitively prove that there is not a population of low 
temperature cores with weak (1,1) lines as well, but we 
do not consider this likely, partly for the theoretical rea- 
son that low density cores would not self-shield against 
radiation and thus would be warmer. 

This does not imply that the central temperatures of 
all these cores are above 9 K. Our model fits a single tem- 
perature to the NH3 spectra. This is a good, but imper- 



fect model. We will overestimate kinetic temperatures 
for cores with a significant temperature gradient due to 
temperature variation along the line of sight as well as 
structure on the sky which is smaller than our beam. 
The contribution from the warm material will tend to 
dominate over the cold. Thus, our results are not incon- 
sistent with cores bei ng colder than 9 K in the cen ter, as 
found in TMC-lc by Schnee & Goodman (20051 where 
the central temperature was 6-7 K in the central 0.03 pc, 
a radius somewhat smaller than our beam. 

Figures [2^ and [2Jd show the temperature distributions 
of our cores broken down into categories and are typical 
of the histograms displayed throughout this paper. Here 
the dashed-blue bars are starless (think empty/cold) 
while the solid-red bars are cores associated with a pro- 
tostar (think full/hot). Certain identifications are dis- 
played in heavy or darker lines, while the lighter blue and 
pink portions correspond to data points with uncertain 
assignment. The two main categories (in this case star- 
less/protostellar) are overlain on each other, while the 
uncertain points are stacked on top of the bars rather 
than being displayed as additional layers. Objects are 
included in the plots if it is possible to derive a kinetic 
temperature from the ammonia spectra, which requires 
a detection of both the (1,1) a nd the (2,2) transition. If 
the (2,2) line is not detected, Rosolowsky et al. (20081 



provides an upper limit on the temperature, 'i'hese are 
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Fig. 1. — Tifin versus the integrated intensity of the (1,1) Une. 
Good solutions for T]jin are shown as plus symbols and 3 a upper 
limits on Tj^in (derived in llosolowsky ct al. 2008; using the noise 
in each individual spectrum) are shown tor the rest of the objects. 



not included on this plot, but see Figure [T] 

Virtually every point below 10 K is starless, but the 
majority of protostellar cores are also relatively cool, be- 
tween 10 and 15 K. They are thus slig htly cooler than 
the 15 K assumed in Enoch et al. (20081. Surprisingly, 3 
of the 4 hottest cores ('i'kin > 20 K) are starless, so not 
all starless cores are cold. These objects are, however, 
at low column density, with relatively weak detections of 
both t he (1,1) and (2,2) line (see Figurejl I. Schnee et al. 
( 2008 1 provide a detailed comparison of these same iNHs 



temperatures with dust temperatures for starless cores in 
Perseus and study the correlation of these temperatures 
with other variables, but are unable to isolate any vari- 
able which strongly influences the temperature of starless 
cores other than that isolated cores are colder than ones 
in clusters. 

The distribution of kinetic temperatures is non- 
Gaussian, though this may be in part because the 
SCUBA/Bolocam surveys which selected most of our 
cores are flux limited, which corresponds to a combined 
temperature/mass limitation. In addition, as noted 
above, we are biased against being able to solve for the 
temperature for very cold objects due to non-detection 
of the (2,2) transition. The isolated cores appear to 
form an approximately normal distribution centered at 
10 K, while the cores in clusters have a roughly uni- 
form distribution between 10 and 15 K. Of the high- 
temperature outliers mentioned above, only two are in 
clusters. There remain two isolated starless cores with 
temperatures above 20 K. 

Kinetic temperature is also a good example of why the 
ANOVA formalism is useful for distinguishing between 
the influence of two categorical variables. The general 
trend is towards colder temperatures in both starless 
cores (Figure [2k) and isolated cores (Figure |2]d) , and we 



might reasonably expect both the presence of a proto- 
star and the surrounding cluster environment to warm a 
core. This variable also illustrates how overlapping these 
populations are: variations within populations are larger 
than the differences between them. 

For the ANOVA test we applied a logarithmic trans- 
formation, which made the distribution closer to nor- 
mal. However, the resulting variable was still long-tailed 
as visualized on a Normal Q-Q plot (not shown). We 
found that cluster environment is highly significant (P 
= 6x10^^), while stellar content is somewhat less signif- 
icant (P=0.01). Clustered cores are 1.16 times warmer 
than isolated ones (95% Confidence Interval:[1.08,1.24], 
while protostellar cores are 1.11 times warmer than star- 
less cores (95% Confidence Interval:[1.01,1.23]). The 
combination of cluster environment and a protostar does 
not increase the temperature beyond the additive in- 
crease from these two individual factors. 

4.1.2. Excitation Temperatures: Tex 
Excitation temperatures (T ex) are only reliably fit a t 



moderate optical depths (see Rosolowsky et al. 20081 



For this sample, the fit temperatures were normally dis- 
tributed with no significant outliers. Stellar content is 
highly significant, while cluster environment is not. Ex- 
citation temperatures are 1.26 K higher in protostellar 
cores (95 % Confidence Interval: [0.73,1.79] P = 6x10"^^). 
Figure [2]3 shows that the protostellar cores have larger 
variance than the starless cores, a fact confirmed by Lev- 
ene's test at the P = 0.02 level. Nonetheless, we remain 
confident that the protostellar cores do have higher ex- 
citation temperatures and propose that this difference 
arises because protostellar cores have higher volume den- 
sities. 

4.2. Volume Densities (n) 
We can make an estimate of the volume density for 



each core u nder the assumption o f LTE from Tc 
and T as ( jSwift et al.||2005[ [CaselU et al.||2002| ) 



T 



(J(rex) - J(rcmb))fcJ(Tkin) 
/i:/l,l(l -J(Tex)) 



X ricrit X P, 



where 



J(T) 



fc(l 



kT 



kin 1 



(3) 



(4) 



and Tcmb — 2.73 K. The escape probability, / 3, is roughly 
estimated for as slab as /3 = (1 — e~'^)/r (Swift et al. 
2005[ ), an estimate which works out to roughly 0.5 for 
our typical values of r. We take t as the maximum r 
(estimated as T(i i) x 0.233 since 0.233 is the maximum 
opacity in any of the hyperfine comp onents). We use 
Ucrit = 2x10^, a fairly standard value (Myers||l999 1. 



A second estimate of the volume density comes from 
simply dividing the mass of a core by the volume, as 
estimated from the radius. This is the approach taken 



in [Enoch et al. ( 2008 ) , though we m odify those masses 
based on our estimates of T^m (see § 4.6 1. We compare 
these densities to those estimated fronTNHa in Figure [3] 
These two estimates systematically disagree, in the sense 
that the Bolocam density estimates are roughly a factor 
of ten larger, and show very little correlation with each 
other. There are plausible explanations, however. First, 
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Fig. 2. — Histograms of kinetic and excitation temperature for ammonia cores broken down by protostellar content and location within 
a cluster. Panels a and b show starless (blue/dashed) and protostellar (red/solid) cores. Uncertain identifications are identified by lighter 
colors and are stacked on top of the certain identifications. Panels c and d show isolated (green/dashed) cores and cores inside of clusters 
(purple/solid). Arrows show median values of the distributions with the thicker arrow pointing to protostellar and cluster means. 



we only see NH3 emission from certain regions of a core, 
so it is not surprising that we see d ensities close to the 
characteristic density (Swade 19891. Second, the Bolo- 
cam masses may be systematically wrong due to uncer- 
tainty in the value of the dust emissivity at 1.1 mm, k. 
Third, the critical density is calculated as 



Aui 

Cul 



(5) 



where the Einstein coefficient fo r the (1,1) transition 



is well established asl.67xl0^' s 
but Gui (the coUisional cross section) is 

known. Calculations of its value range fro 

3 



1 ( |HofcTownes|1983[ ), 



cm' 
s-i 
for 



^ in|Danby et al.| 1 


1988k 


Machin & Kouetf 


(2005 



only poorly 

<io-" 

(for H2) to 5x10^12 cni3 
(for He, both values are 



core-average densities, while the NH3 density calculation 
is for the portion of the core within the beam. As most 
Bolocam cores are larger than our GBT beam and we 



typically point at the center (and presumably densest 
portion) of the core, this will tend to decrease the Bolo- 
cam densities on average. Fifth, Bolocam observations 
partially filter out extended dust emission from core en- 
velopes which is still detected in NH3. Taken together, 
these five factors can easily explain the systematic offset 
and lack of correlation in these two density estimates. 

Regardless of which density determination is more ac- 
curate, protostellar cores appear to have higher volume 
density than starless cores. Nonetheless, because of the 
imprecision of this estimate, we do not attempt further 
statistical tests on this variable. 



4.3. Total Column Density: N(H2) 

We calculate the column density of dust and use this to 
estimate total column density. This is done by measuring 
the flux at the appropriate pixel in the Bolocam map and 
converting to N(H2) using 
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Fig. 3. — Comparison of volume density estimates for protostel- 
lar (red/filled) and starless (blue/open) cores by two independent 
methods. The ordinate shows an estimate from comparing Tj^in 
to Tex from the NH3 data. The abscissa shows an estimate from 
the masses and radii derived from the Bolocam survey of Enoch 
|et al.| | ["2006[ |. The line shows equality. Despite the disagreement, 
both rnethods agree that protostellar cores are denser than starless 
cores. 



N{H2) 



(6) 



where 5']^'=^'" is the flux per beam, ^Ia is the beam sohd 
angle, /ZHa is the mass per hydrogen molecule, is the 
dust opacity, and Bi,{T) is the Planck function. We used 
AiH2 = 2.8, A = 1120 Mm, 6Ihprw = 31" , and = 0.0114 
[cm^ g~^] ( Ossenkopf fc Henning|l994 1 . Th is value of Hy 
corresponds to the wavelength-interpolated [Ossenkopf &:| 
HenningI (|1994 1 value for dust grains with thin ice man- 



tels in dense regions (n = 10^ cm This opacity dif 



fers from that adopted in Andre et al. (19961, which uses 
an opacity a fac tor of 2 lower, consist ent with the rec- 
ommendation of Henning et al. (1995) for lower (< 10^ 
cm~^) density regions. As discussed above, the densi- 
ties of these cores is uncertain, and certainly changes 
throughout the core, making the opacity adopted a sig- 
nificant source of uncertainty. With these values, column 
densities scale like 

The Bolocam map covers virtually all our ammonia 
pointings. Although 10 cores fall outside the Bolocam 
survey, 7 of these are non-detections. For consistency, 
points without a Bolocam flux are not included, even 
if an estimate of the column density is available from 
another source. 

This calculation requires a temperature, for which we 
assume the NH3 kinetic temperature. In general, we ex- 
pect strong coupling between dust and gas temp erature 
at such high densities ( Goldsmith & Langer 1978 1 . In the 
case of protostellar cores, there may be strong tempera- 
ture and chemistry gradients weakening the connection 
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Fig. 4. — Histograms of column density at ammonia pointings as 
derived from the Bolocam flux at the position of the GET point- 
ing and temperatures derived from NH3. Panel a shows starless 
(blue/dashed) and protostellar (red/solid) objects. Uncertain iden- 
tifications are identified by lighter colors and are stacked on top of 
the certain identifications. Panel b shows isolated (green/dashed) 
cores and cores inside of clusters (purple/solid). Arrows show me- 
dian values of the distributions. Gaussian fits are shown and tab- 
ulated in Tabled] 



between ammonia and dust temperatures. Errors in T 
from this gas-dust coupling assumption will propagate 
according to Eqn. [7] 

For the small number (37) of objects where we had a 
NII3 (1,1) detection but were not able to determine a 
temperature from the ammonia (due to an undetected 
(2,2) line), we assume T = 11 K, which is roughly the 
mode and the median of the values in Figure [2j 

The distribution of measured column density is non- 
Gaussian, but becomes essentially normal under the log- 
arithm transformation, reminiscent of the log-normal dis- 
tributi on of column densities throughout the c loud as a 
whole (Ridge et al.'2006| [Goodman et al.||2008| ), but sig- 
nificantly denser. These values generally represent the 
distribution of the highest density peaks in the cloud. 
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TABLE 4 

Gaussian Fits to Logarithmic Column 
Density Distributions 



Category 


Mean Log{N{H2)) 


u Log(N(H2)) 


All 


21.90 


0.32 


Starless 


21.76 


0.34 


Protostellar 


22.07 


0.36 


Isolated 


21.85 


0.35 


Cluster 


22.00 


0.33 



as these are the regions bright enough to be seen in the 
Bolocam map used to select the majority of our targets. 
We fit Gaussians to these distributions in Figure |4] and 
tabulate the fit parameters in Table |4] We use the con- 
version 

N{H2) = 9.4 X 102°cm-2(Av/mag) (8) 



from Bohlin et al. (19781 to convert our N(H2) values 
into an Ay for comparison. 



In the study Goodman et al. (20081, a mean value of 
Ky = 1.18 was found for an extinction map of the whole 
cloud. The mean Ay for our entire sample is 8.4 for 
this fit, confirming that this distribution represents the 
high-density portions of the cloud. 

In our ANOVA formalism we see similar results. We 
derive lower column densities for starless cores than for 
protostellar cores (Figure |4]) , but find no significant dif- 
ference in clustered versus isolated environs. The median 
protostellar core has 2.82 times higher column density 
than the median starless core (95% Confidence Inter- 
val:[1.92,4.16] P = 3xl0~^°). Cluster environment has a 
statistically insignificant influence on dust column den- 
sity. 



4.4. Chemistry 

The derivation of NH3 and CCS co lumn densities are 
explained in Rosolowsky et al. (20081, and obviously re- 
quire detections of the line in question. As the Bolocam 
beam is nearly the same size as the GBT beam (31"), 
we derive a fractional abundance by simply dividing the 
species' column density by N(Il2). The dominant uncer- 
tainty in these chemical calculations is a possible mis- 
match between Tkin and Tdust which would produce er- 
rors in our estimate of N(Il2) of the form described in 



4.4.1. Fractional Abundance of NH3: X-NH3 

We detect NII3 toward 178 of our 193 pointings, so the 
influence of missing the cores with the lowest abundances 
of NII3 (a form of censorship bias) is probably relatively 
small. While the fractional abundance of ammonia does 
not appear to be a strong function of stellar content , clus- 
tered cores generally are less abundant in NH3 (Figure [sji 
& c). Presumably it is the combined radiation environ- 
ment in the cluster that is destroying NH3 rather than 
the influence of an individual embedded protostar. The 
ANOVA test verifies this distinction. Under a logarith- 
mic transformation, the initially skewed fractional abun- 
dance of NH3 remained somewhat skewed. Only cluster 
environment proved to be significant, in the sense that 
the median clustered core had 56% the fractional NII3 
abundance of the median isolated core (95% Confidence 
Interval: [0.35,0.90], P = 0.02). 



4.4.2. Fractional Abundance of CCS: X-CCS 

In the case of CCS, censorship bias becomes a signifi- 
cant problem. Less than half our objects have a reliable 
CCS detection (94 out of 193), and there are no proto- 
stellar cores in clusters with CCS detections (Figure [sj) 
& d). In fact, only five cores within a cluster have a 
CCS detection, and they are all quite weak. In order 
to present a comparable analysis, we assigned Icr upper 
limits to the fractional abundance of CCS based on the 
noise in each individual spectrum, the temperature, and 
an estimate of the width of the line from the linewidth for 
NH3 at that pointing. Then, treating these upper lim- 
its as data points alongside our genuine detections, we 
ran the ANOVA test. The fractional abundance of CCS 
became roughly normal after a logarithmic transforma- 
tion. In this case, both category variables were highly 
significant (still without a significant interaction effect). 
We estimate that the median protostellar cores had 37% 
the fractional CCS abundance of the median starless core 
(95% Confidence Interval: [0.22,0.62] P = 2x10"^), while 
the median clustered core had 32% the fractional CCS 
abundance of the median isolated core (95% Confidence 
Interval: [0.21,0.48] P=8xl0-®). This method overesti- 
mates the CCS abundance in clusters and protostellar 
cores, since more of these objects have 1 a upper lim- 
its (which will typically be higher than the true values). 
Thus, the fractional abundance of CCS in clusters and 
protostellar cores is probably less than presented here, 
and even more distinct from the fractional abundances 
of CCS in isolated and starless cores. 

This distribution is one case in which the division of 
cores into "influenced" and "not-influenced" by our as- 
sociation with 24 micron sources (see Appendix A) pro- 
duces an interestingly different result. The contrast in 
fractional abundance of CCS is more pronounced, with 
virtually all "non-influenced" cores lying at high abun- 
dance. This could be due to radiation from YSO influ- 
encing the chemistry in the outer layers of cores, even if 
that source is not actually embedded within the core. 

4.5. Linewidths 

The total linewidth observed can be decomposed into 
a portion from thermal motion and a portion from non- 
thermal motion. These non-thermal motions may be due 
to chaotic turbulence or ordered bulk motion (e.g. rota- 
tion or shear). In general, all our cores are quiescent, 
in the sense that the non-thermal velocity dispersion is 
less than the thermal component. Nonetheless, it is this 
non-thermal component that is most interesting for un- 
derstanding the dynamics of the cores, and removing the 
thermal component is required for making contrasts be- 
tween different populations, since otherwise the differ- 
ences in temperature would dilute or disguise the results. 

Our non-thermal linewidths are calculated as 



and 



f^thcrmal (NH3) 



f^thcrmal 

(CCS) 



kin 



17m 



H 



kin 



(9) 



(10) 



where 17 and 56 are the molecular weights of NH3 and 
CCS respectively. Our linewidths are measured in terms 
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Fig. 5. — Histograms of fractional abundances of CCS and NH3 broken down by protostcUar content and location within a cluster. Panels 
a and b show starless (blue/dashed) and protostellar (red/solid) cores. Uncertain identifications are identified by lighter colors and are 
stacked on top of the certain identifications. Panels c and d show isolated (green/dashed) cores and cores inside of clusters (purple/solid). 
Arrows show median values of the distributions. 



of the Gaussian cr, so we simply subtract out the thermal 
components for each species to find 



CTnon-thcrmal — ^ total ^thermal' ^^^) 

where crtotai is the measured a of the two different species. 
Our linewidths for NH3 are corrected for optical depth, 
but this was not possible for CCS. 

This calculation also uses the temperature calculated 
for the NH3 line to estimate the non-thermal width of 
CCS. This assumption will be wrong in the case where 
the bulk of the CCS emission comes from a different phys- 
ical region than the NH3. In particular, we might expect 
the CCS emission to emanate from a less dense and more 
turbulent surrounding envelope which is less shielded and 
thus warmer than the very cold and dense region where 
NH3 is found (but CCS is depleted). Indeed, a compar- 
ison of non-thermal linewidths (Figure |6| reveals larger 
non-thermal linewidths in the CCS, suggesting that most 



of the emission from these two molecules arises from 
physically distinct regions. If the CCS comes from a 
warmer portion of the core, we would underestimate the 
thermal component of the CCS line and hence overes- 
timate the non-thermal component. However, due to 
the relatively large mass of CCS, the thermal linewidth 
is small (0.04 km/s) at T = 10 K, and generally much 
smaller than the non-thermal component, so this appar- 
ent difference in non-thermal linewidths is probably due 
to the CCS envelope being genuinely more turbulent. 

4.5.1. Non-thermal linewidth (a"non-thermai^ of NH3 

Figure[7]shows the thermal and non-thermal linewidths 
for NH3 and confirm that most of the cores in our sample 
are quiescent. Figure [8] shows non-thermal linewidths for 
NH3 and CCS divided into the same classes as before. In 
comparing starless to protostellar cores, the non-thermal 
linewidth of NH3 is larger in protostellar cores. The non- 
thermal linewidth for NII3 is normally distributed, and 
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Fig. 6. — Non-thermal linewidths in CCS versus NH3 for starless 
and protostellar cores. Non-thermal motions are larger in CCS 
than NH3, indicating that these molecules are tracing different 
physical regions within the core. Histograms of this plot projected 
onto the axes are presented in Figures [sb and b. 
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Fig. 7. — Thermal linewidths for H2 versus non-thermal 
linewidths leftover in the NH3 profile for starless and protostellar 
cores. Lines demarcate cores with non-thermal linewidths equal 
to thermal and cores with non-thermal linewidths half of thermal. 
Most cores are quiescent; their non-thermal linewidths are less than 
thermal. A histogram of this plot projected onto the ordinate is 
presented in Figure |8^. 

stellar content is significant at the P = 4xlO^^ level. The 
sense is that protostellar cores have more non-thermal 
linewidth (i.e. are more "turbulent") than starless cores, 
though the difference is only 0.025 ± 0.011 km/s, or just 
at the resolution limit of our survey. As turbulent cores 
are less likely to collapse and form stars, the increased 
linewidth in these cores is more likely to be an effect 
rather than a cause of star formation, perhaps indicating 
that feedback from protostars are churning up their natal 



cores. We further illustrate this point in Figure [9] by 
showing the ratio of non-thermal to thermal linewidths 
in NH3, again divided by class. 



Caselli fc Myers] f lMS 
cores within Orion li, Nt 



show that for a sample of 6 
13 linewidth decreases with dis- 
tance from a cluster. Although our clustered cores are 
somewhat more turbulent, we do not see a statistically 
significant difference in non-thermal linewidths between 
clustered and isolated cores. This suggests that the in- 
fluence of relatively weak clusters such as IC348 and 
NGC1333 in Perseus may be significantly weaker than 
the influence of massive clusters such as those in Orion. 

4.5.2. Non-thermal linewidth ('(Tnon-thcrmaiy' of CCS 

For the more diffuse CCS envelope the situation is 
somewhat less clear - starless cores are more prevalent 
at both low (cr < 0.15 km/s) and high (cr > 0.25 km/s) 
linewidths. Note that the x-axis scale is different in the 
two sides of Figure [8j the non-thermal linewidths in CCS 
are generally larger than in NH3. This is also true on 
a case-by-case basis. As illustrated in Figure [8] there 
does not seem to be a significant difference between star- 
less and protostellar cores when it comes to how much 
greater the CCS non-thermal linewidth is than the NH3. 
The non-thermal linewidth of CCS is not normally dis- 
tributed, even after a logarithmic transformation. In any 
event, neither cluster environment nor stellar content ap- 
peared to have a significant influence. No significant 
difference is seen in the clustered/isolated comparison 
(panels c and d) , although the lack of CCS detections in 
clusters precludes any strong statement. 

4.6. Masses and Stability 
The core mass distributions within Perseus have al- 



ready been studie d extensively (Hatchell &: Fuller 2008 
Enoch et al.|2008 l. Our contribution of a better temper 



ature estimate for each core adds only a small correction 
to the previous studies which have had to assume a tem- 
perature for each core (typically lOK for prestellar and 
15K for protostellar). These authors have also distin- 
guished between prestellar and protostellar mass distri- 
butions, but have not looked at the difference between 
clustered and isolated regions. Since we use Tkin from 
NH3 to estimate Tdust in these masses, in cores where 
this assumption does not hold we will have mass errors 
in the same direction and magnitude as for N(H2). See 
94.31 for details. 



Enoch et al. (20081 looks at the stability of prestel- 
lar ari d protostellar cores in Perseus, as does Kirk et al.| 
( 2007 1 using their N2H+ survey in combination with the 



SUUBA cores. We carry out much the same a nalysis, 
but w i th the two following modifications to the |Enoch| 



et al. (20081 results. First, we correct all dust masses 
by using our NH3 temperatures, so our core mas ses are 
related to th ose derived from the Bolocam map in |Enoch| 
et al.|(|2008|) via 



dust 



M, 



,12.85/T 



1 



Enoch 



gl2. 85/10 _ I ' 



(12) 



Second, we have a different division between protostel- 
lar and starless cores, as explained in §2.3[ None of 
these modifications substantially change the conclusions 
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Fig. 8. — Histograms of non-thermal linewidths for CCS and NH3 broken down by protostellar content and location within a cluster. 
Panels a and b show starless (bluc/dashed) and protostellar (red/solid) cores. Uncertain identifications are identified by lighter colors and 
are stacked on top of the certain identifications. Panels c and d show isolated (green/dashed) cores and cores inside of clusters (purple/solid). 
Arrows show median values of the distributions. 



of Enoch et al. ( 2008 ) , which are also broadly in agree- 
ment with the parallel (but independent) dataset s used 
i n the studies of |Hatchell & Fuller (2008) and Kirk et al. 
([2007| 



We adopt the same a ssumption for the ra dial distri- 
bution within a core as Enoch et al. (2008), namely a 



power-law density profile with p ^ r so that we cal- 
culate a = Mdyn/Mdust as defined in Bertoldi & McKee 
(fT992| with 



G ' 



(13) 



where G = 1/232 in units of parsecs, solar masses, and 
km/s. cTdyn is the combination of thermal linewidth for 
H2 and non-thermal linewidth of NH3 added in quadra- 
ture. R is the radius of the Bolocam core as given in 



Enoch et al. (2008); initial indications from a GBT map- 



pmg project of so me of the Bolocam cores in Perseus 
( Pineda et al. 2009 1 suggests that the extent of NH3 emis- 



sion is comparable to the extent of Bolocam emission, 
and thus that the two methods are tracing similar ma- 
terial. This justifies our combination of Bolocam extent 
and mass with NH3 linewidths. 

The virial estimate depends on Tkin in two ways. 
However, since our NH3 spectrum provide high-quality 
linewidths and temperatures, we expect very small er- 
rors on our non-thermal linewidths. Thus the dominant 
dependence on Tkin is through Mjust in the denominator. 

4.6.1. Masses: Mdust 

Figure [lO] shows the distribution of core masses for 
both starless/protostellar and clustered/isolated. The 
distribution is roughly log-normal, though likely infiu- 
enced on the low-mass side by incompleteness; the 50% 
completeness limit for the average sou r ce siz e (70") is 
at 0.8 M0 as estimated in Enoch et al. (20061. We find 
that cluster environment is insignificant, but protostellar 
cores have 1.32 times the median mass of starless cores 
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Fig. 9. — Histograms of the ratio of non-thermal to thermal 
linewidths for NH3 broken down by protostellar content and lo- 
cation within a cluster. Panel a shows starless (blue/dashed) and 
protostellar (red/solid) cores. Uncertain identifications are identi- 
fied by lighter colors and are stacked on top of the certain identi- 
fications. Panel b shows isolated (green/dashed) cores and cores 
inside of clusters (purple/solid). Arrows show median values of the 
distributions. 



Fig. 10. — Histograms of core masses broken down by protostel- 
lar content and location within a cluster. Panel a shows starless 
(blue/dashed) and protostellar (red/solid) cores. Uncertain iden- 
tifications are light colors and are stacked on top of the certain 
identifications. Panel b shows isolated (green/dashed) cores and 
cores inside of clusters (purple/solid). Arrows show median values. 
The 50% completeness limit is shown for the average source size 
(70"). 



(95% Confidence Int erval:[l. 00,1.76] P = 0.05 ). This is 
consistent with what Hatchell & Fuller ( 2008[ ) found for 
cores in Perseus using the SCU15A data, where a K-S test 
found a 4% chance that the masses of protostellar and 
starless cores were drawn from the same distribution (as- 
suming lOK for starless and 15K for protostellar). Since 
NH3 temp eratures for protostellar c ores are typically less 



than 15 K, Hatchell & Fuller (20081 are underestimating 



the masses for protostellar cores relative to ours, draw- 
ing the two populations somewhat closer together. This 
explains why our test, despite having more degrees of 
freedom, provides roughly the same level of confidence 
that these populations are different. Levene's test re- 
veals that the variances of the starless and protostellar 
distributions are unequal at the P = 0.01 level. 



4.6.2. Vmal Parameter a: Mdyn/Mdust 

Figure[TT]shows the distribution of the virial parameter 
(a) for both starless/protostellar and clustered/isolated. 
Again, this variable requires a logarithmic transform to 
make it normal, and our test shows no significant influ- 
ence from either stellar content or environment. Most 
of our objects are gravitationally bound (a < 2) and 
roughly a third are in gravitation al virial eq uilibrium 
(a < 1). This is in agreement with |Enoch et al., ( |2008"] ) 
(which is a similar an alysis on the sam e datasets), but 
in stark contrast with iKirk et al. ( 2007 1 , who find via a 
similar analysis of N2H"*" and SCUBA cores that no cores 
in Perseus are in gravitational virial equilibrium and few 
are gravitationally bound. 

Rather than attempt a detailed reconciliation between 
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these two studies, which rely on different algorithms to 
define clumps in maps with different noise and sensitivity 
limits, we cross-match our Bolocam sources with SCUBA 
positions within 60"and compare average properties of 
the two populations of cores. This method allows us to 
u nderstand why our cores are all bound while the cores 



Kirk et al. (20071 are all unbound, but unfortunately 
does not answer the question of which method better 
determines whether a core is truly bound. The largest 
difference comes from the masses and sizes of the cores 
defined from the thermal emission maps; the non-thermal 
component of the N2H^ lines are very similar to the 
NH3 non-thermal linewidth s. Compared to o ur Bolocam 
cores, the SCUBA cores of Kirk et al. (2007) are: much 
less massive (62%) and somewhat smaller (24%), with 
slightly higher total (temperature-adjusted) linewidths 
(110%). Since 

a^{Ra^)/M (14) 

if these SCUBA cores were scaled by these median values, 
a would change by a factor of 2.7 and the results would 
be in rough accord with our own - mostly bound cores 
with a significant number virialized. 

5. COMPARISON WITH OTHER SIMILAR 
SURVEYS 

Apart from the other studies carried out in Perseus 
which we discussed above, w e also compare our results 
to those of Jijina et al. ( 1999 ), long the standard compila- 
tion of NH3 core properties, and the NH3 and extinction 
study of the Pipe Nebula w hich was constru cted in a 
fairly similar fashion to ours ( Lada et al.||2008 ). 

We prov ide a compar i son w ith just one of the re- 



sults from Jijina et al. ( 1999 1 which illustrates some 



of the adv antages our uniform survey provides. iJl 



jina et al. ( [1999) find that cores within clusters are 
more massive (as estimated from their NH3 lines assum- 
ing virial equilibrium) and have a higher ratio of non- 
thermal to thermal linewidths. By comparison, all our 
objects fit w ithin a narr ow r ange in Figure 4 from | Jijina 



et al. (1999J (see Figure 12), and there does not appear 
to b e significant distinc tions between the sub-samples. 



The Jijina et al. 

log (m; 



( 1999 1 sub-sample medians run from 

= 0.7 to 1.5, and O'non-thermal/o'thermal 

0.6 to 1.6. As shown in Figure [12] 

log 



runs from from 

all our objects are clustered near the lowest mass 
{Myir /Mq) between 0.1 and 0.4) and least turbulent ob- 
jects (t Tnon-thermal / "thermal between 0.5 and 0.8). 

The jj^ ina et al.l (|1999|) sample includes cores from 
Perseus, most of which (15/20) are considered to be pro- 
tostellar and isolated. Our objects within Perseus have 
similar levels of turbul ence, but a re sig nificantly less mas- 

T wo sim ple ex- 
(1999) assume 



sive than those of j Jijina et al. ([i999|) 



planations for this are that Jijina et al 



a distance of 350 pc for Perseus, where we use 250 pc 
and that the observations of Perseus in this other sam- 
ple were taken largely from Effelsberg (40"beam) and 
Haystack (88"beam), the latter of which in particular 
is much larger than our 31"beam, and thus might eas- 
ily blend two sources which we resolve. Both these ef- 
fects would tend to increase the masses. The much larger 
range in core prop erties identified in the full sample of 



Jijina et al. ( 1999 1 could be the result of cores being in- 
tluenced by larger scale (i.e. cloud) properties. 
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Fig. 11. — Histograms of the virial parameter, a. Panel a shows 
starless (blue/dashed) and protostellar (red/solid) cores. Uncer- 
tain identifications are lighter colors and are stacked on top of the 
certain identifications. Panel b shows isolated (green/dashed) cores 
and cores inside of clusters (purple/solid). Dashed and dotted lines 
show the criteria for gravitationally bound and in gravitational 
virial equilibrium. Vertical arrows show median values. 



Lada et al. 



^ (|2008| conclude that for the Pipe Nebula, 

a nearby molecular ^loud with little star formation, most 
cores are thermally-dominated objects bound only by ex- 
ternal pressure. This contrasts with our cores, which ap- 
pear to be gravitationally bound. The Pipe core masses 
come from an extinction map, and their dynamical in- 
formation comes mainly from a pointed survey of C^^O 
(objects were examined for NH3 but with a low success 
rate). Due to the lack of NH3 detections 
(2008) assumed 10 K for every core 



Lada et al. 



It is possible that the objects detected in the Pipe are 
fundamentally d ifferent objects than the millimeter cores 
studied herein. [Alves et al.j (2007) note that the cores 
identified within thcH^ipc have typical densities between 
5x10^^ cm^"^ and 2x10^ cm^^^, much lower than the den- 
sities of the Bolocam cores discussed herein. The low 



16 



Foster et al. 



1.5 - 



o 



0.5 - 



0.0 



-I— 1 — I — I — I — I — r 



] — I — I — I I I I I I — I — I — I — I — I — I I I I I I I 



IRAS ond Cluster 



1 — I — I — r- 



Perseus Doto 
(j;jino,1999) 




No IRAS 
but Cluster 



IRAS but no Cluster 



No IRAS ond 
no Cluster 



f- 



Perseus Doto (this study) 
I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ 



0.4 



0.3 



s' 0.2 



o 



0.1 



0.0 



— 1 1 1 1 1 1 1 1- 

Perseus Doto (this study) 



No Stor 
but Cluster 



No Stor 
No Cluster 



Stor ond Cluster 



Stor but no Cluster 



0.0 0.5 



1.0 
o 



1.5 2.0 

non - thermol/^ thftrmol 



2.5 3.0 



0.4 



0.6 

<^ non- 



0.8 



non-thermol/ ^thftrmol 



1.0 



Fi g. 12. — Median an d quartile log (M„ jr/M(ri) vs. c Trtr^rt—t hpi-m!,] /o'thetTmi,] for the four different sub-samples of objects as in Figure 4 
from [jijina et aL] | |1999l l. On the left, the |Jijina et al.] ( |1999^ results (IRAS point sources were used as a proxy for protostellar), which 
spa n a much larger ra nge in both axes and exhibit more separation between classifications than our Perseus data. Cores in Perseus from 
the |Jijina et al. | | |1999[ l database are typically more massive. We zoom into our data on the right, illustrating that our different classes of 
objects largely overlap in this diagram. 



rate o f NH3 detections in the study of [Rathborne et aT] 
( 2008 1 appears to confirm that these cores are relatively 
low density. The Pipe cores could be either younger cores 
(which will eventually become denser) or a distinctly dif- 
ferent population of objects which would be undetected 
by millimeter surveys at the distance of Perseus. In the 
latter scenario, the cores we define in Perseus represent 
j ust a portion of the true core mass function. 

Lada et al. (2008) do a least-squares fit of Log(Mass) 
versus Log(a) an d find a slope of -2/3, in ag reement with 
the prediction of Bertoldi & McKee ( 1992 1 for pressure- 
confined cores. \Ve perform a similar experiment on our 
sample of obj ects. Our masses s pan a similar range, but 
Lada et al.l (12008!), almost all of our cores 

This 



in contrast to 

are consistent with being gravitationally bound 



could be due to genuine differences between the two re- 
gions under consideration as discussed above or due to 
the different linewidths of the tracers used. Linewidths 
are the obvious source of systematic differences simply 
because of dominance of systematic uncertainty of this 
term {a ^ cr^yn)- C^'^O traces a lower density envelope 
than NH3, and thus using it to estimate the linewidths 
of cores will produce an overestimate of the non-thermal 
component. This overestimate was on average a factor 
of 4 within the Pipe cores observed in both tracers. 

An unweighted least-squares fit to our full data set 
provides the same slope (M ^ a~^^^^), however a closer 
look reveals some cautions with this analysis. First, 
we calculate the uncertainty on our derived parameters, 
propagating the er rors on masses and radii quoted in 
Enoch et af] ( |2006t and the errors on tem peratures and 



^ (|2008| into our de- 

i'his ignores systematic errors in the 



linewidths from Kosolowsky et al. 
rived quantities, 
mass estimate. With these estimates we are able to state 
that the low-mass cores are almost all consistent with 
being gravitationally bound without external pressure. 



and that high-mass cores are unstable to gravitational 
collapse (unless supported by magnetic fields). 

To consider the influence of selection effects, we take 
the completeness limits on Bolocam cores from ' Enoch] 
[Ct al._ (2006j . One constraint is simply due to the beam 
size, which prohibits detection of objects with FWHM 
< 34". There are also curves in mass and radius de- 
rived from Monte Carlo experiments on simulated data 
and simulated reduction routines which arc roughly flux 
limits (high mass objects are difficult to detect if too ex- 
tended) - we use the 90% completeness limit presented. 
These curves are transformed into a and mass constraints 
which requires an estimate of (Tnon-thcrmai- Our ability 
to measure the linewidth is generally not a constraint 
within our survey, but since our completeness in alpha 
is thus dependent on three variables (mass, radius, and 
linewidth), it is not possible to represent it faithfully in 
this plot. 

To provide a useful guide, we use simple upper and 
lower bounds for linewidths. The flux limit provides 
an upper bound on alpha (we can't detect things too 
extended/far-from-bound), so we use the 75% quartile 
of the non-thermal linewidths for the protostellar (0.174 
km/s) and starless (0.148 km/s) sub-samples when cal- 
culating the total linewidths of an object. This provides 
a reasonably conservative number upper bound on the 
linewidths and thus alpha. The beam-width limit is a 
lower bound on alpha (we can't measure a radius less 
than the beam), so we use the 25% quartile linewidths 
for the full population. Figure [T3| shows the data points, 
errors, and completeness estimates for starless and proto- 
stellar cores. Our completeness estimates nearly bound 
the data points on the upper side, and could contribute 
to the shape of this relation. 

In particular, completeness effects could be severely 
limiting our ability to detect self-gravitating cores at 



Properties of Dense Cores 



17 



masses below 1 Mq. At high-masses, the constraints 
are not as severe. That is, we can be reasonably confi- 
dent that there are no high-mass cores which are grav- 
itationally unbound, and the highest mass cores are all 
estimated to be significantly self-gravitating (i.e. likely 
undergoing collapse) . 

6. CONCLUSION 

As summarized in Table [3] the median values of cores 
in clusters are statistically significantly different from the 
median values of isolated cores in the following ways: 

1. Higher kinetic temperatures (Tkin) 
Cluster Cores = 12.9 K 
Isolated Cores = 10.8 K 

2. Lower fractional abundances of CCS 
Cluster Cores = 0.6 xlQ-^ 
Isolated Cores = 2.0x10^9 

3. Lower fractional abundances of NH3 
Cluster Cores = 1.2x10-^ 
Isolated Cores = 2.9x10^^ 

The median values of cores associated with protostars 
are statistically significantly different from the median 
values of starless cores in the following ways: 

1. Higher kinetic temperatures (Tkin) 
Protostellar = 11.9 K 

Starless = 10.6 K 

2. Higher excitation temperatures (T^x) 
Protostellar = 7.4 K 

Starless = 6.1 K 

3. Higher column density (N(H2)) 
Protostellar = 1.2x10^^ cm'^ 
Starless — 0.6x10^^ cm^^ 

4. Lower fractional abundance of CCS 
Protostellar = 1.4x10"^ 

Starless = 2.4 xlO'^ 

5. Higher non-thermal motion (cTnon-thermai NH3) 
Protostellar = 0.14 km/s 

Starless = 0.11 km/s 

6. Higher masses 
Protostellar =1.5 Mq 
Starless = 1.0 M© 

Neither cluster environment nor protostellar content 
makes a significant difference to the dynamical state of 
cores as estimated by the virial parameter. Most cores 
are bound. 

We are not able to assess whether differences in core 
properties arise from the influence of protostars/clusters 
or whether protostars and clusters form around cores 
with a given set of properties. Nonetheless, we propose 
explanations for these trends which present a coherent 
picture. 

We see evidence that the cluster environment has a 
significant influence on the ~ 10"*^ cm~^ gas characteristic 
of NH3. This is most straightforwardly interpreted as 
the influence of radiation, which could (1) increase the 
temperature and modify the chemistry, decreasing the 



abundance of both (2) CCS and (3) NH3. T his sort of 
behav ior is seen, for example, in the models of |Lee et al.| 
( [2004| in which CS (which behaves like CCS) and NH3 
abundances are suppressed in the outer regions of a core 
by the presence of a UV field. The slightly higher column 
density seen in clustered cores could be responsible for 
lowering the CCS abundance (via depletion) , but can not 
easily explain the lowered abundance of NH3 . 

The small temperature increase seen is not sufficient to 
be distinguishable in the virial analysis we have carried 
out. We see no evidence that the mass distribution of 
cores is different between clusters and isolated regions, al- 
though further study of this question in warranted. This 
suggests that low mass clusters may form most of their 
stars in a mode fundamentally similar to isolated star 
formation. Other forms of feedback, such as outflows or 
shells may be important in determining the properties 
of cores in clusters, but we are unable to point to any 
specific evidence of such an influence. 

Clusters with higher-mass star formation may present 
a better chance to study the effects of feedback and ex- 
amine the differences between low- and high- mass star 
formation. Although the relatively low-mass s tars in 
NGC 1333 are lik ely destroying the clouds (see Guter- 
muth et al.||2008| and references therein), these mid-H 



stars still have a relatively weak impact on their envi- 
rons compared to the 0-stars which dominate a region 
such as Orion. Furthermore, the clusters in Perseus have 
between 100 and 300 stars, while th e Orion Nebula Clus- 
ter has ^ 1500 (Porras et al. 112003 1. Studies on the am- 



monia clumps in Orion already suggest that NH3 cores 
here have different properties from these seen in Perseus 



(Har iu et al.||1993|]Caselli fc Myers]|1995 | 'Wilson et al 
yyyl. Ihe cores m Orion from Jijina et al. (1999) show 



a large range in many variables, with some clear distinc- 
tions (Tkin is greater in clusters than in isolated cores), 
but larger samples would be useful. A study which con- 
tained a large number of NH3 cores both near and far 
from the clusters in Orion would provide definitive state- 
ments about the influence of massive clusters on dense 
cores. 

We propose that some of the properties of protostellar 
cores are due to the presence of a protostar, while other 
properties are indicative of the environment necessary for 
star formation. Higher non-thermal linewidths and tem- 
peratures would seem to inhibit rather than promote star 
formation. Since we see both these features in protostel- 
lar cores we propose that we are seeing the influence of 
the protostar on its natal core. Chemical changes would 
seem likely to be either caused by star formation or due 
to an independent variable such as time (since it takes 
time for cores to both evolve chemically and to collapse 
into a star). 

On the other hand, it has been suggested that pro- 
tostars only form in cor e s with sufficiently high column 
density Johnstone et al. (20041. We interpret our higher 



excitation temperatures as arising from higher physical 
densities, which might promote star formation (for in- 
stance by changing the thermal physics as in [Larson 
2005|). Und er the competitive accretion model ( [Bonnell 



et al.|200ip , cores which are already going to form a star 



might be expected to draw in additional mass, but it is 
also plausible that more massive cores are simply more 
likely to collapse. Thus we are left with the conclusion 




Fig. 13. — Log masses and virial parameter a for our sample of cores broken down into starless (blue/open diamonds) and protostellar 
(red/stars). The horizontal broken black lines show the conditions for gravitationally bound (dotted) and virial equilibrium (dashed). 
Colored lines show estimated completeness limits using the 75% quartile linewidths for these two samples at the 90% (dot-dashed) flux- 
detection limit. The solid black line shows an estimate of the influence of beam size in detecting compact sources. Completeness could be 
significantly influencing this relation. 



that some of the differences in protostellar cores are due 
to the influence of those young stars, while other differ- 
ences may explain why those particular cores are chosen 
out of the ensemble to birth a star. Our lack of distinc- 
tion between starless and protostellar cores on the ques- 
tion of dynamical stabihty seems to suggest that most 
of our starless cores are relatively long-lived rather than 
transient structures. 

Ultimately, a theory of star formation should be able 
to reproduce the properties of the cores presented here, 
both pre- and post- the star formation event, both in- 
side clusters and in isolated regions. Matching all these 
properties will provide much tighter constraints on the 
physics of star formation than simply matching the mass 
function of these cores. Much progress has already been 
made by numerical simulations of star formation within 
a molecular cloud. To match the observations herein. 
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7. APPENDIX A: CLASSIFICATION OF CORES AS STARLESS OR PROTOSTELLAR 

In principle, the question we want to ask when classifying a dense core as protostellar or starless is: "Is a star 
which formed within this dense core still embedded within the core?" Evolution presents an intractable complication 
to answering this question. A dense core may survive the formation of a star. This star is then subject to different 
forces (it does not feel gas pressure, for instance) and may move away from the core. At some point we will no longer 
be able to associate this YSO with its natal core. This problem may be particularly intense in clusters where we must 
also worry about another problem, that a neighboring star may pass close to a different dense core and thus appear to 
be associated. Of course, if we wish to measure the generic impact of a young star upon a core, associating this star 
with this new core may still be acceptable. 

Even putting aside that problem, dense cores form part of a smooth larger distribution of material and thus may 
have many legitimate different definitions of their extent . There are thus different valid definitions of what constitutes 
an embedded YSO. This is complicated by the fact that our analysis depends on two different data sets, a millimeter 
continuum map for which a measure of the 2-dimensional extent may be made, and pointed NH3 observations, for 
which no such measure can be established. Indeed, the situation is worse, for some of our NH3 pointings are at 
positions for which no sensible extent can be measured from the millimeter map, either because the pointing is outside 
the Bolocam map, or because the millimeter continuum at a position has too low a signal-to-noise ratio. To maintain 
our sample size and limit additional biases we seek a general metric for whether the NH3 gas observed at a given 
position is associated with recent star formation or not. 

We combine the results of two previous studies which have made this distinction with two of our own methods. The 
t wo literature s t udies have more sophisticated methods, but rely on data not available for our entire survey. 



Enoch et al. ( 2008 1 provides a classification of Bolocam cores as starless or protostellar based on constructing a 
catalog of co l d VsO s. This catalog is a subset of the candidate young stellar objects (YSOc) from the c2d catalog 



(Evans et al. (20071) with some additional very red objects included. Then, using the sizes for Bolocam cores defined 



Enoch et al. (2006 1, they consider a core to be protostellar if there is a cold YSO within 1 FWHM of the center of 



the core. This method provides a reasonably clean definition for millimeter cores, but it is not extensible to objects for 



which we are unable to define a millimeter core. Furthermore, some of the cores in Enoch et al. (2006| are quite large 
(mean minor-axis FWHM is 58"± 17" and mean major-axis FWHM is 80"± 27", with major axes~FWHM extending 
out to 150"). Our GBT beam is probing gas within about 30"of the center of these cores, so while a YSO may 
legitimately be impacting the total flux from a large millimeter core, it may not have any impact on the dense gas we 
a re tracing. 



Hatchell et al. (20071 also builds a catalog of young YSOs (class or 1) by considering SCUBA-identified cores 
(supplementing this list with additional Bolocam cores) and constructing and fitting SEDs from these long-wavelength 
fluxes with many additional shorter wavelengths including 2MASS and IRAC fluxes. The association of shorter 
wavelength fluxes with the (sub)millimeter poin ts and the inclusion o f short wavelength (2MASS) points in constructing 
the SED are two essential differences with the Enoch et al. ( 2008 1 results. The two methods agree reasonably well, 
but disagree on 10 out of the 76 cores which both classify. We consider every core where these methods disagree to be 

an uncertain classification. 

Our own two methods are rougher but applica ble to all sources. First, we take the the cold YSO list of 'Enoch et al.| 
( [2OO8 ) and the YSOc list from the c2d catalog (Evans et al. (2007)) and at each GBT pointing we see whether one of 



these objects falls within 1 FWHM of the beam (30"). i'his method shows good agreement with the two classifications 
from the literature, disagreeing with one of these methods only 13 times (and 10 of these times the two literature 
methods disagree, as noted above) 



Sec ond, we examine the proximity and strength of 24 /im sources from the c2d MIPS data for Perseus (RebuU et al. 
20071. Some 24 /im sources are background galaxies, but the majority are YSOs. By determining the average spatial 
density of 24/im sources as a function of intensity we are able to make an estimate for each position: "Is there a 24/im 
source close and bright enough that it is unlikely to be there by chance?" Not including the ambiguous cases discussed 
above, this method categorizes 43 more cores as being protostellar than the other methods, largely due to the fact 
that cores are very clustered, which places many starless cores in close proximity to a bright 24 /im source which 
nonetheless lives in another core. There is also significant bright, diffuse 24 /im emission towards the east of Perseus. 

The 24 /im method was very useful in assessing positions where neither literature method provides a classification. 
Since these are necessarily weaker cores (or outside the area mapped), it is not surprising that they are overwhelmingly 
starless. A final case of note is the mini-cluster L1455, which lies on the edge of the Bolocam and c2d surveys. Here 
we were restricted to the 24 /im criteria and other literature information about the YSOs in this well-studied region. 
There were 4 cores considered in this fashion. We did use the uniform 24 /im criteria to divide cores into two categories 
"influenced" and "not-influenced" and produced the plots within this paper under this division, though we do not show 
these plots. In all cases the same qualitative behavior was observed, in the sense that starless cores are mostly "not- 
influenced" and protostellar cores are mostly "influenced" . We mention the one ca se wh ere this distinction produces 
substantially different distributions within the text (fractional abundance of CCS, [4.4.21. 

In summary, if any two of the three main methods disagree we examine the images my hand and consider the 
reasons for the disagreement. This is normally due to a dispute over what the extent of the core really is, or whether 
a somewhat distant YSO should be considered associated, but sometimes hinges on the quality of the SED used to 
make the classification. We make an assessment and assign this object a classification, but mark it as uncertain. These 
objects are all presented with a brief comment in Table |5] 
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TABLE 5 

Classification of objects as Protostellar/Starless 







Bolocam 


Class 


Hjnocn 


xlatcnell 


CZa YbLIC 


24 jxm 


Comments 


# 


detected? 


core? 








in beam? 




1 


N 


N 


L 






N 


N 


Outside of Bolocam and SCUBA. 


2 


Y 


N 


L 


N 




N 


N 




3 


Y 


Y 


L 


N 




N 


N 




4 


Y 


Y 


L 


N 




N 


N 




5 


Y 


Y 


L 


N 




N 


Y 




6 


Y 


Y 


L 


N 




N 


N 


Outside of SCUBA. 


7 


Y 


Y 


s 


Y 


Y 


Y 


Y 




8 


Y 


N 


S? 




Y 


N 


Y 


Hatcliell lists a YSO, but close to NH3SRC 7. 


9 


Y 


Y 


L 


N 




N 


N 


Outside of SCUBA. 


10 


Y 


N 


L 






N 


Y 




11 


Y 


Y 


L 


N 




N 


N 


Outside of SCUBA. 


12 


Y 


Y 


s 


Y 


Y 


Y 


Y 




13 


Y 


Y 


L 


N 




N 


N 


Outside of SCUBA. 


14 


Y 


Y 


s 


Y 


Y 


Y 


Y 




15 


Y 


Y 


L 


N 


N 


N 


Y 




16 


Y 


Y 


L 


N 




N 


N 


Outside of SCUBA. 


17 


Y 


Y 


L 


N 


N 


N 


N 




18 


Y 


Y 


L 


N 




N 


N 




19 


Y 


Y 


L 


N 




N 


N 




20 


Y 


Y 


L 


N 




N 


N 




21 


Y 


Y 


L 


N 




N 


N 




22 


Y 


Y 


s 


Y 


Y 


Y 


Y 




23 


N 


N 


L 






N 


N 


Outside of Bolocam and SCUBA. 


24 


Y 


Y 


L 


N 




N 


N 




25 


Y 


N 


L 






N 


N 


Outside of SCUBA. 


26 


Y 


N 


L 






N 


N 


Outside of SCUBA 


27 


Y 


N 


L 






N 


N 


Outside of Bolocam and SCUBA. 


28 


Y 


N 


L 






N 


N 


Outside of Bolocam and SCUBA. 


29 


Y 


Y 


L 


N 




N 


N 




30 


Y 


N 


L 






N 


N 




31 


Y 


Y 


S 


Y 


Y 


Y 


Y 




32 


Y 


Y 


S 


Y 


Y 


Y 


Y 




33 


Y 


N 


L 




N 


N 


Y 




34 


Y 


Y 


S 


Y 


Y 


Y 


Y 




35 


Y 


Y 


S 


Y 


Y 


Y 


Y 




36 


Y 


N 


L 






N 


Y 


At edge of SCUBA and Bolocam. Starless in 


















literature.' 24 /im source in NH3SRC37. 


37 


Y 


N 


s 






Y 


Y 


38 


Y 


N 


L 






N 


N 


Outside of Bolocam and SCUBA. Starless in 
literature.' 


39 


N 


N 


L 






N 


N 




40 


Y 


Y 


s 


Y 


Y 


Y 


Y 




41 


Y 


Y 


L 


N 


N 


N 


N 




42 


Y 


Y 


S? 


Y 


N 


N 


Y 


Hatchell has a starless core slightly to the north 
but Enoch and c2d have a source on edge of GBX 
beam. 


43 


Y 


Y 


s 


Y 


Y 


Y 


Y 




44 


Y 


Y 


s 


Y 


Y 


Y 


Y 




45 


N 


N 


L 




N 


N 


N 




46 


Y 


Y 


s 


Y 


Y 


Y 


Y 




47 


Y 


N 


L? 




N 


Y 


Y 


c2d source just within beam. 


48 


Y 


Y 


S 


Y 


Y 


Y 


Y 




49 


Y 


Y 


L 


N 




N 


Y 




50 


Y 


Y 


S? 


N 


Y 


N 


Y 


Small Bolocam core with Enoch just outside 
beam. Hatchell has YSO. 


51 


Y 


Y 


L 


N 




N 


Y 




52 


Y 


Y 


s 


Y 




Y 


Y 




53 


Y 


Y 


s? 


Y 




N 


Y 


Big Bolocam core. Nearby c2d source outside 
beam. 


54 


Y 


N 


L 






N 


N 




55 


Y 


N 


L 






N 


N 




56 


Y 


Y 


s 


Y 




Y 


Y 


Multiple c2d sources within beam and Bolocam 
core. 


57 


Y 


N 


L 






N 


Y 




58 


Y 


Y 


S 


Y 


Y 


Y 


Y 




59 


Y 


Y 


s 


Y 


Y 


Y 


Y 




60 


Y 


N 


L 






N 


N 


Outside of Bolocam and SCUBA. 


61 


Y 


N 


L? 






N 


Y 


Enoch object just outside beam. 


62 


Y 


N 


L 






N 


Y 




63 


Y 


N 


L 






N 


N 




64 


Y 


Y 


S 


Y 


Y 


Y 


Y 




65 


Y 


Y 


s 


Y 


Y 


Y 


Y 




66 


Y 


Y 


s 


Y 


Y 


Y 


Y 




67 


Y 


Y 


s 


Y 


Y 


Y 


Y 
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TABLE 5 — Continued 



NH3SRC'' 


NH3 


Bolocam 


Class*' 


Enoch'^ 


Hatchell'^ 


c2d YSOc"^ 


24 /im 


Comments 


Tr 


detected? 


core? 








in beam? 




68 


Y 


N 


S 




Y 


Y 


Y 




69 


Y 


Y 


L 


N 


N 


N 


Y 




70 


Y 


N 


L? 




Y 


N 


Y 


Nearby Enoch sources are inside other cores. 
These fluxes probably contribute to Hatchell 

classification. 


71 


Y 


Y 


S 


Y 


Y 


Y 


Y 




72 


Y 


Y 


L? 


N 


Y 


N 


Y 


Hatchell source (#62) has an oddly rising SED. 
No c2d source within beam. 


73 


Y 


Y 


L? 


Y 


N 


N 


Y 


Large Bolocam core but all c2d sources are out- 
side of beam and inside other cores. 


74 


N 


N 


L 






N 


N 




75 


Y 


Y 


S 


Y 


Y 


Y 


Y 




76 


Y 


Y 


S 


Y 


Y 


Y 


Y 




77 


Y 


Y 


s 


Y 


Y 


Y 


Y 




78 


Y 


N 


s 




Y 


Y 


Y 




79 


Y 


Y 


L? 


N 


Y 


N 


Y 


Hatchell source (#70) has flat/rising IRAC 
fluxes. Not too odd, but no c2d source within 
beam. 


80 


Y 


Y 


L 


N 


N 


N 


Y 




81 


Y 


Y 


S 


Y 


Y 


Y 


Y 




82 


Y 


Y 


s? 


Y 


N 


Y 


Y 


Closest Enoch source in neighboring core, but c2d 
source at center of beam. 


83 


Y 


Y 


L 


N 


N 


N 


N 




84 


Y 


Y 


S 


Y 


Y 


Y 


Y 




85 


Y 


N 


L 






N 


Y 




86 


Y 


Y 


L 


N 


N 


N 


N 




87 


Y 


Y 


S 


Y 


Y 


Y 


Y 




88 


Y 


Y 


L 


N 


N 


N 


Y 




89 


Y 


Y 


S 


Y 




Y 


Y 


Outside of SCUBA. 


90 


Y 


N 


L 






N 


N 


Outside of SCUBA. 


91 


Y 


Y 


S 


Y 


Y 


Y 


Y 




92 


N 


N 


L 






N 


N 




93 


Y 


Y 


L 


N 




N 


Y 




94 


N 


N 


L 






N 


N 




95 


Y 


Y 


S 


Y 


Y 


Y 


Y 




96 


Y 


Y 


L 


N 




N 


N 




97 


Y 


Y 


L 


N 


N 


N 


N 




98 


Y 


N 


L 






N 


N 




99 


Y 


Y 


S 


Y 


Y 


Y 


Y 




100 


N 


N 


L 






N 


N 


Outside of Bolocam and SCUBA. 


101 


Y 


N 


L 






N 


N 




102 


Y 


N 


L 






N 


N 




103 


Y 


Y 


S 


Y 


Y 


Y 


Y 




104 


Y 


Y 


L 


N 


N 


N 


N 




105 


Y 


Y 


S 


Y 




Y 


Y 




106 


N 


N 


L 






N 


N 




107 


Y 


Y 


L 


N 




N 


N 




108 


Y 


Y 


L 


N 


N 


N 


Y 




109 


Y 


Y 


L 


N 




N 


N 




110 


Y 


N 


L 






N 


N 




111 


Y 


Y 


L 


N 




N 


N 




112 


Y 


Y 


L 


N 




N 


N 




113 


Y 


Y 


L 


N 


N 


N 


Y 




114 


Y 


Y 


L 


N 




N 


Y 




115 


Y 


N 


L 






N 


Y 




116 


Y 


Y 


L 


N 




N 


N 




117 


Y 


Y 


S 


Y 




Y 


Y 




118 


Y 


Y 


S? 


Y 


N 


Y 


Y 


Hatchell says starless (# 82) with IRAC upper 
limits, but Enoch and c2d have a centered source. 


119 


Y 


Y 


s 


Y 


Y 


Y 


Y 


Two c2d sources. 


120 


N 


N 


L 






N 


N 




121 


Y 


Y 


s 


Y 


Y 


Y 


Y 




122 


Y 


N 


L 






N 


N 




123 


Y 


Y 


S 


Y 


Y 


Y 


Y 




124 


Y 


Y 


L 


N 




N 


Y 




125 


Y 


Y 


L 


N 




N 


N 




126 


Y 


Y 


S 


Y 


Y 


Y 


Y 




127 


Y 


Y 


L 


N 




N 


N 




128 


Y 


Y 


S? 


N 




Y 


Y 


Starless in Enoch, but close/faint 24/im source 
and c2d source. 


129 


Y 


N 


L 






N 


N 




1.30 


N 


N 


L 






N 


N 




131 


N 


N 


L 






N 


N 




132 


Y 


Y 


L 


N 




N 


N 


Outside of SCUBA. 


133 


N 


N 


L 






N 


N 
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1\TTJ 
INM3 


rJoiocam 




tnocn 


rlatcneli 


CZQ loUC 


24 fiv[i 


Oomments 


M- 

W 


detected? 


core? 








in beam? 






134 


N 


N 


L 






N 


N 




135 


N 


N 


L 






N 


N 




136 


N 


N 


L 






N 


N 




137 


N 


N 


L 






N 


N 




138 


N 


N 


L 






N 


N 




139 


Y 


N 


L 






N 


N 


Outside of SCUBA. 


140 


N 


N 


L 






N 


N 




141 


Y 


Y 


L 


N 




N 


N 


Outside of SCUBA. 


142 


Y 


Y 


L 


N 


N 


N 


Y 


Nearby 24/im source. 


143 


N 


Y 


S 


Y 


Y 


Y 


Y 




144 


Y 


Y 


L 


N 




N 


N 




145 


Y 


Y 


L 


N 


N 


N 


Y 


Diffuse 24^m 


146 


Y 


Y 


L 


N 




N 


Y 


Diffuse 24^m 


147 


Y 


Y 


L 


N 


N 


N 


Y 


Diffuse 24/xm 


148 


Y 


N 


L 






N 


N 




149 


N 


N 


S 






Y 


Y 


c2d source right at center. 


150 


Y 


Y 


L 


N 




N 


Y 




151 


Y 


Y 


L 


N 




N 


Y 




152 


Y 


Y 


L 


N 




N 


Y 


Unclear 24/im detection. 


153 


Y 


Y 


L 


N 




N 


Y 


Unclear 24/zm detection. 


154 


Y 


Y 


L 


N 




N 


Y 


Unclear 24/im detection. 


155 


N 


N 


L 






N 


Y 


Diffuse 24fim 


156 


Y 


Y 


L 


N 


N 


N 


Y 




157 


Y 


Y 


L 


N 


N 


N 


Y 




158 


Y 


Y 


L 


N 




N 


Y 




159 


Y 


Y 


L 


N 




N 


Y 




160 


Y 


Y 


S 


Y 


Y 


Y 


Y 




161 


Y 


Y 


S 


Y 


Y 


Y 


Y 




162 


Y 


Y 


s 


Y 


Y 


Y 


Y 




163 


Y 


Y 


L 


N 


N 


N 


Y 


Unclear 24/im detection. 


164 


Y 


Y 


s? 


Y 


N 


Y 


Y 


Enoch has a source in NH3SRC 164-165-169 clus- 
ter and two other c2d sources within this beam. 


165 


Y 


Y 


L? 


N 


N 


Y 


Y 


Enoch source is within 164, although at edge of 
beam. 


166 


Y 


Y 


L? 


Y 




N 


N 


Enoch says protostellar, but others methods say 
starless, and Bolocam core is very diffuse. 


167 


Y 


N 


L 






N 


N 




168 


Y 


Y 


S 


Y 




Y 


Y 




169 


Y 


Y 


L? 


Y 


N 


N 


Y 


This Enoch source is in NH3SRC 164. 


170 


Y 


Y 


L 


N 


N 


N 


N 




171 


Y 


Y 


L? 


N 


Y 


N 


N 


Hatchell says protostellar, but everything else 
says no. 


172 


N 


N 


L 






N 


Y 




173 


Y 


Y 


s 


Y 




Y 


Y 




174 


Y 


Y 


L 


N 


N 


N 


Y 




175 


N 


N 


L 






N 


Y 




176 


Y 


Y 


L 


N 


N 


N 


Y 




177 


N 


N 


S 






Y 


Y 


Two c2d sources within this core. 


178 


Y 


Y 


s 


Y 


Y 


Y 


Y 




179 


N 


N 


L 






N 


Y 




180 


Y 


Y 


L 


N 


N 


N 


Y 




181 


Y 


Y 


L 


N 




N 


N 




182 


Y 


N 


L 






N 


N 




183 


Y 


Y 


L 


N 


N 


N 


Y 




184 


N 


Y 


S 


Y 




Y 


Y 


c2d source at center. 


185 


N 


N 


L 






N 


N 


Outside of SCUBA. 


186 


N 


N 


L 






N 


Y 


Outside of SCUBA. Pointing misses star in tail 
of B5 


187 


N 


N 


L 






N 


N 


Outside of SCUBA. 


188 


Y 


Y 


L 


N 


N 


N 


Y 




189 


Y 


Y 


L? 


Y 


N 


N 


Y 


Hatchell defines a starless core. Larger Bolocam 
core includes NH3SRC 192. 


190 


Y 


N 


L 






N 


N 




191 


Y 


N 


L 






N 


Y 




192 


Y 


Y 


S 


Y 


Y 


Y 


Y 




193 


N 


N 


L 






N 


N 


Outside of SCUBA. 



These reference numbers are from 
- Protostellar, L? — Uncertain Starless, S? 



Rosolowsky et al. 



(2008 



Uncertain Protostellar 



which c an be consulted for in dividual core position s and propertie s 

' " 'd' ' ' " 



Enoch et al. 



I 2008 I 



Hatchell et al, 
T 



associated if a YSO candidate falls within 1 FWHM of the GBT beam (3l") of the pointing location 



' 20071 



Evans et al. 



Jijina et al. 



I 1999 I 



■ Starless, S 



\ 2007 I, considered 
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8. APPENDIX B: SENSITIVITY TO RANDOM AND SYSTEMATIC TEMPERATURE ERRORS 

Many of the distributions we measure and describe depend critically on a temperature derived from modeling the 
NH3 spectra in order to convert observed quantities into physical units. This dependence on temperature can be 
divided into two cases: 

1) The distributions which use only the NH3 and CCS data (Tox, Tkin, and both non-thermal linewidths) require 
that our determination of temperature accurately measures the temperature of the bulk (H2) gas. The most serious 
error possible is a systematic one. 

2) The distributions which also rely on the Bolocam maps (N(H2), Mdust, Q^, X(CCS) and X(NH3)) require that T^in 
from NH3 adequately approximates Tdust • This is an estimate and may be biased in some way. We show in Appendix 
B that the best estimate for this discrepancy has a small influence on our results. Nonetheless, for this half of the 
analysis the probabilities quoted are valid only under the assumption that T^in for NH3 adequately approximates 

Tdust ■ 

As the two simplest quantities to assess, we perform experiments on the non-thermal linewidth of NH3 (case 1) and 
the total column density (case 2). We perform 4 distinct Mo nte Carlo experiments: 



a) We apply the random uncertainties of T^jn derived in Rosolowsky et al. (2008) to both case 1 and case 2 for 
10,000 individual realizations and calculate the percentage of cases in which our distmctions hold. The test we perform 
on the result is to simply take the median of the ratios for the protostellar and starless populations, discarding the 
cluster information for simplicity. Furthermore, the confidence intervals and P- values derived from ANOVA are not 
calculated using measurement errors (instead, they look at the width of the sample distributions). Thus, we do not 
expect our simulations to fully reproduce the confidence intervals present in ANOVA. As shown in Figure |14| the 
addition of these random errors has a negligible influence on our results. 

b) We apply our best estimate of the systematic offset between T^in and Td^st to case 2 (1 K underestimate of 
the dust t emperature) and again run 10,000 realizations with random uncertainties. We take this estimate of 1 K 
from Galli et al. (2002 1 Figure 3, which shows the dust and gas temperatures within a starless core of the appropriate 

This figure is for the the low radiation case, which we believe to 

As shown in Figure [15] this sort of 



density range to be traced with Bolocam and NH3 

be appropriate for the relatively low intensity star-formation present in Perseus 
systematic offset has a small influence on our result. 

c) We apply a pessimistic estimate of the systematic offset between Tkin and Tdust to case 2. Here, NH3 over- 
estimates the temperature of starless cores by 2 K and under-estimates the temperature of protostellar cores by 2 K 
This will increase the temperature differen ces between starless and pr otostellar cores, approximating the temperatures 



difference assumed in Enoch et al. (2008^ and 



uncertainties. As shown in J^'igure To 



(20071. We again run 10,000 realizations with random 



Hatchell et al. 

sucn an offset severely diminishes (but does not eliminate) the difference between 
column densities towards protostellar and starless cores. 

d) We adopt a systematic error such that Tkin in warm sources is underestimated by an arbitrarily chosen amount 
(New Tkin = (l+t^) Tkin)- We allow S for vary from to 1 and examine at what value of 5 the difference in these 
distributions is eliminated. As shown in Figure |17| significant bias is required to eliminate the difference in column 
density, and the difference in non-thermal linewidths is robust against an error of this sort. 

These tests show that our results are very robust against the influence of random errors on the temperature. This 
is mostly due to the high-quality of the spectra, which allows a very precise estimate of the temperature. The non- 
thermal linewidth result is robust against systematic errors in the temperature determination. Systematic errors in 
the temperature determination can be constructed which eliminate the distinctions we see in quantities which rely on 
the dust emission map, particularly if these systematic differences work in different directions for the two populations 
(i.e. protostellar/starless). However, under our best estimate of the possible systematic error (Trial b) our results are 
affected by only a small amount. 
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1.8 1.9 2.0 2.1 2.2 

Medlon Protostellor/Storless N(H2) 
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a 
E 

3 



2000 - 



1000 - 





1.15 




1.16 1.17 1.18 1.19 1.20 1.21 1,22 
Medlon Protostellor/Storless <7„„„.„^,^| 



Fig. 14. — Results of 10,000 simulations adding estimated random uncertainties to Tj^in for N(H2) (left panel) and for the non-thermal 
NH3 linewidth (right panel). In both cases the median ratio of the protostellar and starless populations are shown. On the left, we see 
asymmetric scatter around the original value of 2, but no point approaches 1 (which would be the case where our statement that protostellar 
cores are at higher column density would be incorrect). On the right, a much smaller relative variation is present in non-thermal linewidths 
as this variable is less sensitive to changes in Tjjjn. The original value of this ratio is 1.184. 
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c 
o 
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1000 - 



800 - 



o 
£ 



0) 

E 

3 



600 - 



400 - 



200 - 




K7 
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2,2 



Fig. 15. — Results of 10,000 simulations adding estimated random uncertainties to Ti^;,, for N(H2) and a 1 K underestimate of the dust 
temperature. The median ratio of the protostellar and starless populations is shown. The bias shifts the ratio from the original ratio 
of 2 down to peak near 1.9 with some asymmetric scatter. No point approaches 1 (which would be the case where our statement that 
protostellar cores are at higher column density would be incorrect). 
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o 
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£ 



0) 

E 

3 



1500 - 
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1.15 



Fig. 16. — Results of 10,000 simulations adding estimated random uncertainties to Tjjin for N(H2) and the assumption that NH3 
overestimates the dust temperature of starless cores by 2 K and underestimate the dust temperature by 2 K in protostellar cores. The 
median ratio of the protostellar and starless populations is shown. The original ratio of 2 has been shifted substantially, and this distribution 
now approaches 1. For larger systematic errors of this particular kind, the distinction between the populations is erased or reversed. 
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40 60 
Percent bios in T 



100 




40 60 
Percent bios in T 



100 



Fig. 17. — The ratio of N(H2) and non-thermal lincwidth for our two populations under the assumption of a systematic problem in using 
NH3 to determine tiic bulk gas temperature of the form Tjjin = (1+5) T^kin- In both cases the median ratio of the protostellar and starless 
populations are shown. For N(H2) this systematic error needs to reach 60% before significantly influencing the ratio, and the ratio does 
not reach equality until 80%. The result for non-thermal linewidths are almost entirely robust to this sort of error. 



